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I. INTRODUCTION 


1.1 Motivation . 

The purpose of this study is to examine the theory of sound pro- 
pagation in a viscous, heat-conducting fluid, initially at rest and in 
a uniform state, and contained in a rigid, impermeable duct with iso- 
thermal walls. 

In the past ten years, research in duct acoustics has been largely 
motivated by the necessity of reducing the noise level associated with 
the commercial operations of modern jet aircrafts. A large portion of 
this noise is internally generated in the jet engine itself, and is 
primarily due to the rotating turbomachinery blades and the combustion 
process. This noise then propagates in the form of acoustic modes 
through the inlet and exhaust pipes. Even though, there has been a 
great deal of effort devoted to the understanding of basic aerodynamic 
noise generation mechanisms, one of the most efficient ways of minimi- 
zing the radiated sound still consists of absorbing as much of the 
acoustic energy as possible along the propagation path inside the inlet 
or exhaust pipe. This is essentially achieved by treating the duct 
walls with a suitable sound absorbing material. However, the viscous 
and heat-conducting properties of the medium itself also contribute to 
the attenuation of the radiated sound. Most aeronautical applications 
fall in the high-frequency wide-tube range, i.e., the frequencies of in- 
terest and the cross-sectional dimensions of the duct are such that dis- 
sipation due to friction and heat conduction is restricted to a thin 
acoustic boundary layer close to the duct walls, and only constitutes a 
small part of the total attenuation. It is therefore not surprising 
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that most of the current studies in duct acoustics of jet engines assume 
the medium to be inviscid, Nevertheless, acoustic boundary layer atten- 
uation has to be taken into account in any systematic evaluation of the 
total attenuation rate. 

These effects also play a primary role in the propagation of blood 
pressure pulses in human arteries and in the performance of fluid trans- 
mission lines. The medium is then a liquid and heat conductivity is 
usually found to be negligible. The dimensions of the tube and the char- 
acteristics of the fluid are such that, in contrast to the previous in- 
stance, viscous forces are dominant throughout the cross-section, and 
viscous dissipation is large. This is the so-called low -frequency -nar- 
row-tube range and the medium may then no longer be considered inviscid. 

From a theoretical point of view, this problem may be viewed as ''ne 
of the few instances where the basic equations pertaining to the unsteady 
motion of a viscous, heat-conducting, and compressible fluid are amenable 
to analytical treatment. Even thouj^i these equations are linearized, one 
expects the solution to retain most of the essential features of more 
complicated flow situations. The fundamental concern of this study is 
to develop a systematic methodology that would be of use in analyzing 
complex aerodynamic noise problems. 

1.2 Review of the Literature . 

In 1868, Kirchhoff published his famous study, "On the Influence of 
Heat Conduction on Sound Propagation in Gases." At about the same time, 
Regnault (1868) made the first measurements of sound attenuation in the 
sewers of Paris. A detailed account of Kirchhoff 's theory is given in 
Lord Rayleigh's Theory of Sound (1877). Kirchhoff considered both vis- 
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cosity and heat conduction, and derived a general dispersion relation 
pertaining to axisymmetrlc waves in cylindrical tubes. This relation 
was then solved for plane waves In the high-frequency and low-frequency 
approximations. Many papers on this same subject have been written 
since, and it would be tedious to discuss them here in any detail. One 
may, however, within the extensive literature available, distinguish a 
few main trends. 

A large number of investigators have sought to extend Kirchhoff's 
plane mode results to the entire frequency domain. Weston (1953a), in 
particular, computed additional terms in the first-order expansions 
given by Kirchhoff, in both the high- and low-frequency cases, and re- 
vealed the existence of an additional very-high-frequency range. Iber- 
all (1950), Brown (1962), and Rott (1969), neglecting radial pressure 
gradients in the governing equations, derived a solution for the plane 
mode valid in the low and high frequency domains. Finally, Shields et 
al. (1965), and very recently Tijdeman (1969, 1975) made a comparison 
of the previous results with a numerical solution of the dispersion re- 
lation. It is worth mentioning at this point that another group of re- 
search workers, Sexl (1930), Womersley (1954), Lance (1955), and Uchida 
(1956), concerned more specifically with the propagation of sound in non 
heat-conducting liquids, studied oscillating flows in tubes filled with 
an incompressible viscous fluid. Even though there is no wave propaga- 
tion in this case, the axial velocity profile bears a close resemblance 
to the corresponding compressible flow result. In particular, in the 
high-frequency limit, the amplitude of the axial velocity presents a 
maximum close to the duct wall. This phenomenon, experimentally dis- 
covered by Richardson and Tyler (1929) , is usually referred to as the 


Richardson Annular Effect. All theoretical results pertaining to the 
pil'd or fundamental mode were abundantly checked experimentally by Fay 
(1940), Kemp and Nolle (1953), Meyer and Guth (1953), and Weston (1953b). 

From inviscid acoustics, however, one would expect the existence at 
high frequencies, of higher-order propagating modes. A few studies have, 
Indeed, been devoted to the determination of their attenuation character- 
istics. The first attempt seems to have been made by Hartlg and Lambert 
(1950), whose theoretical formulation later proved to be unsatisfactory. 
Cremer (1948) suggested simulating viscous and thermal effects close to 
a rigid wall by an equivalent impedance, and Beatty (1950a, b) and Morse 
and Ingard (1968) then used the well-known soft-walled -duct model to com- 
pute the attenuation rate of higher order modes. Shaw (1950, 1953) di- 
rectly perturbed Kirchhoff*s dispersion relation to recover essentially 
the same results. Correpsouding expressions were also derived for a 
purely viscous fluid by Bogert Q950) , Elco and Hughes (1962), and Cohen 
and Tu (1962). Nayfeh (1973) recently extended Cremer's equivalent im- 
pedance concept to situations where the medium is inhomogeneous and non- 
uniformly moving. Paradoxically, the full problem has not yet been 
solved numerically, although the case of waves in a viscous fluid con- 
tained in a cylindrical tube was treated by Gerlach and Parker (1967) 
and very recently by Scarton and Rouleau (1973). Scarton and Rouleau, 
in particular, used the method of eigen-valleys to show the existence of 
a previously unknown family of vorticity-dominated modes. Their work 
seems so far to be the most comprehensive study of the viscous effects 
on sound propagation in ducts. 

Even though wave propagation in tubes is the ultimate concern of 
the present investigation, it is a particular case of the more general 
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theory cl the smell unsteady motions of a viscous, heat-conducting and 
compressible fluid, and a significant portion of this work has been de- 
voted to a new mathematical formulation of such a class of problems. It 
Is therefore appropriate to review here some of the studies made in this 
wider context. Lagerstrom »t el. (1949) examined the fundamental prin- 
ciples underlying the linearized system of equations governing the flow 
of a viscous compressible fluid, they showed the existence of two types 
of waves, namely, longitudinal or pressure waves, and transversal or 
vorticity waves. In particular, a proof was given for the statement 
that any linearized flow may be split into a longitudinal and a trans- 
versal component. Wu (1956) generalized their results to the case of a 
viscous and heat-conducting fluid. Chu and Kovasznay (1958) were pri- 
marily interested in the measurement of fluctuating quantities in super- 
sonic turbulent flows, and Introduced the concept of three main modes of 
fluctuations, namely, sound, vorticity, and entropy modes. This idea 
was recently expanded in the wider context of aerodynamic noise theory 
by Doak (1973). Finally, two basic papers of general interest and deal- 
ing mainly with sound propagation in an unbounded medium are worth men- 
tioning here. Truesdell (1953) conducted an exhaustive review and cri- 
tique of Kirchhoff's theory of free space sound absorption, for fluids 
with different heat-conducting properties, and for different frequency 
regimes. Lighthill's survey (1956) on finite -amplitude sound waves shed 
further light on the physical mechanisms responsible for sound absorp- 
tion. Both studies have interesting discussions on the controversial 
issue of bulk viscosity. 

A century has passed since Kirchhoff's investigation and it is ap- 
parent from this brief review of the literature that a complete treat - 
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ment of the problem is yet unavailable. In the next section, we pro- 
ceed to a general discussion of the scope and goals of the present work. 

1.3 Scope of This Study . 

This study is essentially concerned with the theory of small per- 
turbations of a viscous, heat-conducting, and compressible fluid. As 
such, the second chapter is devoted to an alternative mathematical for- 
mulation of the problem, suggested by Doak (1973) in the general frame- 
work of aerodynamic noise theory. In this approach, the fluctuating 
velocity field is considered as the superposition of acoustic (related 
to pressure), thermal (related to entropy), and viscous (related to 
vorticity) parts. Such a decomposition reduces the linearized Navier- 
Stokes equations to a system of three partial differential equations 
for three basic unknown function.*: an acoustic scalar potential; a 

thermal scalar potential; and, e vii ous vector potential. No restric- 
tive assumptions are made regarding the nature of the fluid or the mag- 
nitude of the Prandtl number. The physical implications of this formu- 
lation are discussed in detail, as well as particular cases where sig- 
nificant simplifications may be achieved. 

The previous methodology is applied in the third chapter to a study 
of small -amplitude fluctuating motions in a duct of constant width with 
rigid and isothermal walls. The governing system of equations is pro- 
perly non-dimensionalized in terms of five non-dimensional parameters 
which are functions of the cross-sectiona 1 dimension of the tube, the 
frequency of the perturbation and the characteristic properties of the 
fluid. The remainder of the investigation further assumes the duct to 
be two-dimensional, in order to minimize analytical complications. It 
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is then shown that each frequency may be associated with four distinct 
types of wave motion, namely, antiplane or inplane waves of symmetric or 
antisymmetric shapes with respect to the duct axis. Each family is char 
acterized by a specific dispersion relation between a given circular fre 
quenc and the corresponding complex wave numbers. The relations per- 
taining to antiplane waves may readily be solved for arbitrary values of 
the parameters whereas the relations characterizing inplane motions are 
transcendental in nature and require the implementation either of a nu- 
merical scheme or of some approximation procedure. It is found that per 
turbation methods give relatively simple mathematical expressions in 
terms of the characteristic parameters of the problem, and, at the same 
time provide a convenient framework for the interpretation of the main 
physical phenomena occurring during the wave motion. 

Consequently, Chapter IV is concerned with a perturbation study of 
the inplane dispersion relations in three ranges of frequencies or equi- 
valently three ranges of duct widths: the low-frequency-narrow-tube 

range; the high-frequency-wide-tube range; and, the very-high-frequency- 
very -wide -tube range. It is shown that, in addition to the usual pres- 
sure-dominated modes encountered in inviscid propagation problems, two 
other families of entropy- and vorticity-domina ted modes have to be con- 
sidered. Expansions for the attenuation rates, phase velocities, and 
mode shapes are then obtained for each family in each frequency regime. 
The physical implications of the results are then discussed and a quali- 
tative comparison is made with the numerical solution of Scarton and 
i\o uleau (15/3), in the limit of zero heat-conduction. 

In the concluding remarks, we summarize the essential contributions 
of the research and suggest opportunities for further study based on the 

experience gained in the present undertaking. 
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II. THEORETICAL FORMULATION OF THE SMALL 
AMPLITUDE FLUCTUATING MOTIONS OF A VISCOUS, HEAT-CONDUCTING 
AND COMPRESSIBLE FLUID 


2,1 Introduction . 

In problems concerned with the acoustic motions of an inviscid 
fluid initially at rest and in a uniform state, it is customary to de- 
rive from the basic linearized equations of Fluid Mechanics a single 
wave equation for a perturbation velocity potential $(r,t) : 


D 2 <t = 0 


( 2 . 1 ) 


where 



( 2 . 2 ) 


is the wave operator and a^ is the isentropic speed of sound. The 
other dependent variables such as the velocity perturbation V r , and 
the pressure perturbation p* , are then directly related to the po- 
tential $ through the classical expressions: 


V* = grad $ (2.3) 



where is the density of the medium. 

This procedure presents the important advantage of reducing a wide 
variety of acoustic problems to the determination of one scalar field 
$(r,t) obeying the wave equation (2.1) with appropriate initial and 
boundary conditions. 

In this chapter, we wish to generalize this formulation to 
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situations where the medium is viscous and heat-conducting. We will be 
led to introduce, in addition to the acoustic potential defined above, 
two other fields, namely, a thermal scalar potential and a viscous vec- 
tor potential. Governing equations for these "auxiliary" variables will 
be derived, and the physical implications will be discussed in light of 
the mathematical formulation. 

2.2 Linearized Navier-Stokes Equations and Underlying Assumptions . 

Let us consider the small amplitude motions of a simple thermo- 
dynamic fluid initially at rest and in a uniform state. Each dependent 
variable such as pressure, density, velocity, ... may be written in the 
form: 

Q(r,t) = Q q +Q' (r,t) (2 * 5) 

where Q q is independent of r and t and characterizes the initial 
state of the fluid, and Q'(r,t) designates a fluctuating quantity. In 
particular, since V Q = 0 

V (r , t) = V'(r,t) (2.6) 

Such a decomposition is then substituted for each dependent variable 
into the equations of mass, motion, entropy, and the equations of state, 
pertaining to a Newtonian fluid obeying Fourier’s law of heat conduction. 
When all non-linear terms in the perturbation variables are neglected, 
the following linearized equations result: 

Equation of Mass ^ H-p^divV* = 0 (2.7) 

Equation of Motion Pq - grad p ' - curl w'+TJ^grad div V' (2.8) 
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Equation of Entropy 

0 T ~ = k V^T 1 

p o T o at V 1 

(2.9) 

Equations of State 

p' = a 0 2p’ + [|f] S' 

(2.10) 


V s o 



T' - [g] I p' +[§] S' 

“'si , ss p „ , 
p o ,s o p o ,s o 

(2.11) 


In the above relations, Uq is the coefficient of shear viscosity, k Q 
the coefficient of heat conduction, the isentropic speed of sound, 
and 77^ is the dilatational viscosity given by 


% = K 0 + 3^0 


( 2 . 12 ) 


Kq being the coefficient of bulk viscosity. All other symbols have 
their usual significance. 

The system (2.7) -(2.11) constitutes the starting point of this in- 
vestigation and we discuss below the main assumptions which have been 
introduced in the process of its derivation: 

(a) Continuum Hypothesis : The characteristic length scales of the 

particular problem at hand are assumed to be much larger than the mean- 
free-path, so that the fluid may be treated as a continuum. In the case 
of propagation in ducts, the wavelength and duct diameter must both be 
large in comparison with the mean-free-path. 

(b) Homogeneous Medium at Rest : The steady flow variables are in- 
dependent of position and time and the velocity of the medium is identi- 
cally zero. 

(c) The fluid is assumed to satisfy the Navier -Stokes relation be - 
tween stress and rate of strain as well as Fourier's law of heat con - 
duction. 
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(d) Small Amplitude Motion : In other words, fluctuating and steady 


quantities satisfy the following relations: 


-♦ 

11 

« 1 ; 

£l 

« 1 ; 

ill 

a o 


p 0 


s o| 


(2.13) 


This assumption justifies the linearization of the governing equations 
and drastically simplifies the mathematical formulation of acoustic 
problems. 

(e) Laminar Motion . The presence of turbulent motion would at 
once invalidate the linearization of the basic equations since turbu- 
lence is characterized by strong non-linear interactions. Unfortunately, 
very few studies on the transition of unsteady laminar flows have been 
made. However, in the case of oscillatory flow in circular ducts, Ser- 
geev (1966) and Nerem et al. (1972) determined experimentally a criter- 
ion for transition of the form 


& 


.e - const, 
c 


when 


>4 


(2.14) 


In the above relation, Re^ is the critical value of the Reynolds num- 

A 

ber based on the peak velocity U , and k and R^ are non -dimens ion- 
alized parameters which will be introduced in Chapter 3. They are de- 
fined as follows: 


A 

Re - HI 


A 


V, 


0 


j a.d 

k * ~ ; R 

a 0 V V 0 


(2.15) 


where d is the duct diameter, u) the circular frequency of the lami- 
nar flow oscillations, and the kinematic viscosity. 


It will be shown in Chapter IV, that for values of yj kR^ larger 
than four, the flow is characterized by a thin acoustic boundary layer 


close to the duct wall of thickness 


6 



(2.16) 


A crude justification of relation (2.14) may then be given by applying 
the critical Reynolds number criterion for a steady flat plate boundary 
layer to the unsteady boundary layer of thickness, 6 , i.e., 



- 1000 


c 

Upon combining this relation with the estimate of 
the dependence of the critical Reynolds number on 


(2.17) 


6 given by (2.16) , 



may be shown 


to be expressed by (2.14). 

In this study, it is assumed that the peak Reynolds number is lower 
than the critical Reynolds number given by (2.14) so that the unsteady 
flow is laminar. 

(f) Simple Thermodynamic Fluid : Nonequilibrium effects such as 

vibrational relaxation are assumed to be negligible and the fluid is in 
local thermodynamic equilibrium. The thermodynamic state of the fluid 
at any position and time is completely described by two state variables 
only, for instance, density and entropy. A»1 other state variables such 
as pressure and temperature may then be expressed in terms of density 
and entropy, as in the two linearized equations of state (2.10) and 
(2.11). The partial derivatives in these equations may conveniently be 
expressed in terms of five characteristic properties of the medium, 
namely, the isentropic speed of sound a^ , the temperature T^ , the 
density p n , the specific heat at constant pressure c p , and the 
ratio of specific heats . Equations (2.10) and (2.11) may then be 


12 


written in the following manner: 


>>' ■ *o 2 p,+p o*o 


/< y 1)T o 


T' 




(y £^!» p ,*Vo 

C C 


( 2 . 10 )' 


( 2 . 11 )' 


The form of the coefficients in the above relations is derived in Appen> 
dix A, by making use of the first and second law of thermodynamics. No 
additional assumption regarding the nature of the fluid needs to be in- 
troduced. In the case of a perfect gas, for instance, a Q could be re- 
placed by its expression as a function of , but this restriction 
does not have to be made in order to derive (2.10)' and (2.11)'. 

The linearized equations (2.7) - (2.9) and (2.10)', (2.11)' consti- 
tute a system of five linear partial differential equations for five 

— ♦ 

unknown functions p 1 , p* , S' , T* f and V* . We now proceed to the 
definition of the acoustic, thermal, and viscous potentials, and to the 
derivation of a corresponding system of three partial differential equa- 
tions for these three potentials. 


2.3 Acoustic, Thermal, and Viscous Potentials. Governing Equations : 
According to Helmoltz* theorem, the velocity field V 1 can always 
be written as the sum of an irrotational part grad $ , and a solenoidal 

-*4 

part curl A so that 

V* = grad $ 4-curl A (2.18) 

We will further assume that the so-called viscous vector potential A 
satisfies the additional condition 
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div A ■ 0 


(2.19) 


From (2.10)', the density p' may be expressed as a function of p' 
and S' . When (2,18) and the resulting equation for p' are substi- 
tuted Into the equation of mass (2.7), one obtains 



( 2 . 20 ) 


In order to satisfy (2.20) identically, Doak (1973) suggested in his, 
"Momentum Potential Description of Unsteady Fluid Flows," to further 
decompose the scalar potertial $ into two parts: an acoustic poten- 

tial $ giving rise to pressure fluctuations through the relation: 

A 



-Vo 272 *. 


( 2 . 21 ) 


and a thermal potential $ giving rise to entropy fluctuations 
through the relation: 


3S' 

at 




( 2 . 22 ) 


The velocity field V' is now the sum of three distinct acoustic, ther- 
mal, and vortical components, and is of the form: 


V' = grad 4* + grad $ , *!• curl A 
a th 


(2.23) 


From (2.10)' and (2.11)', the other dependent variables may immediately 
be related to $ , $ ^ * an ^ A in the following manner: 





(2.24) 
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°V<V 1 )c p, 


V ^th^V 1 ^ 3 


(2.25) 


Equations (2.21) - (2.25) provide the necessary relations between the 

"auxiliary" dependent variables 4^ , 4 , and A , and the "physical" 

variables, pressure, entropy, velocity, density, and temperature. In 

the derivation of these relations, use has been made of the equation of 

continuity and the equation of state. We now seek to obtain the partial 

differential equations governing 4 , 4 , and A , by requiring that 

s cn 

the remaining equations of motion and entropy be satisfied. Substitu- 
tion of (2.23) and (2.22) into (2.8) and (2.9) respectively, leads to 
the following system: 


gr.d[A(* + 4. th )-^’ 2 (* +* th > + $ 


(2.26) 


+ curl [|^ - v Q V 2 A ] * 0 


V 2 [4 




th Vr 


- T'] 


(2.27) 


where P is the Prandtl number of the medium: 
r 


r k o 


(2.28) 


It is very tempting to replace System (I) by the following simpler 
system of three partial differential equations: 


^-(4 +4 J-~ v 2 (4 +4 . )+ E -^ = 0 

dt v a th' p Q v a th p Q 


(2.29) 


ft ■ -o’ 21 * 0 


'(v n -l)c t 


th a P 
o r 


T' = 0 


(2.30) 


(2.31) 
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Each solution of (II) is obviously a solution of (I). But there are 
solutions of System (I) which are not solutions of System (II). In 
Appendix B, we prove that, to each solution {4^ , 4^ , A*] of (I), 

one may associate a solution {4^ , 4^ , A] of (II), which will give 
the same values for the physical variables V', p' , and S' . It may 
be concluded that Systems (I) and (II) yield the same number of possible 
physical states. 

Equations (2.29) and (2.31) may immediately be expressed solely in 
terms of 4^ and 4 ^ by differentiating them once with respect to 
time and making use of relations (2.21) and (2.25). The final results 
of this section are then formulated in the following "Splitting Theorem." 

Splitting Theorem : If { V ' , p' , S' , p' , T ' } satisfy the li- 

nearized equations of Fluid Mechanics (2.7) -(2.9) and (2.10)' - (2.11)', 
then there exists an acoustic potential 4^ , a thermal potential 4 ^ , 
and a viscous vector potential A such that the physical variables are 
represented by: 


V* « grad $ + grad $ , + curl A 
a th 


_ p a 2 7 2 4 

5t H 0 0 a 


with 



th 


th 


94 


th 


(2.32) 

(2.33) 

(2.34) 

(2.35) 


(2.36) 

(2.37) 
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and where ^(r,t) , $ th (r,t) ^nd A(r,t) satisfy the following system: 


a o u 

(2.38) 

bt P r v th °0 L) P r V a 

(2.39) 

|-‘' 0 7 2 A-0 

(2.40) 


In equation (2.36), equation (2.39) was used to derive an alternative 
3T * 2 

expression for in terras of 4^ only. □ S in equation (2.38) 

is the modified wave operator. 


□ 


2 

S 





p o a o 



(2.41) 


The traditional wave operator as defined in (2.2), is Immediately re- 
cognized in the first two terms of (2.41). The last third-order term 
represents ettenuation brought about by the dilatational viscosity ^ . 
The subscript S has been added to notify that the operator pertains to 
the laentropic speed of sound a^ . 

Furthermore, cross differentiation of (2.38) and (2.39) leads to 
the additional fourth order partial differential equations: 


where 



(2.42) 


(2.43) 


The above operator is a modified wave operator similar to 
in (2.41), but pertaining to the isothermal speed of sound 


0^ defined 

a Q^\fro * 
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Discussio n: Upon examination of the previous theorem, it may be 

noted that the propagation problem has been reduced, as expected, to the 
determination of an acoustic scalar potential ♦ , a thermal scalar po- 

tential ♦ . , and a viscous vector potential A . 

tn 

The acoustic component gives rise to pressure, density, and velocity 
perturbations. It is inherently coupled to the thermal component, as 
seen from (2.38) and (2.39), and therefore is bound to generate entropy 
fluctuations. The thermal component gives rise to entropy fluctua- 
tions and, as a result of coupling vlll generate the acoustic 

component. Both components satisfy equation (2.42) independently. 
However, coupling is still achieved through the lower-order rela-.ons 
(2.38) and (2.39). Finally, the velocity fluctuations associated with 
and ♦ ^ are irrotational so that investigators have often consi- 
dered these two components as one single longitudinal part. The uncoup- 
ling of the vortical component represented by A is the most signifi- 
cant feature of the splitting theorem. The vector potential is gov- 
erned by a diffusion equation end only gives rise to rotational velocity 
fluctuations. However, in the presence of solid surfaces, one expects 
the vortical component to generate the other two through the no-slip 
boundary conditions. 

We wish to emphasize that the decomposition presented here is not 
the only possible one, and several alternative schemes have indeed been 
proposed. In particular, Wu (1956) distinguished a longitudinal and a 
transversal component, without any further splitting of the longitudinal 
part. Chu and Kovasznay (1958) defined three modes, namely, pressure, 
entropy, and vorticity modes, but the first two modes do not coincide 
with the acoustic and thermal components considered in this study. In 
both of these investigations, the governing equations were expressed in 
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terms of the physical veriebles. We cen see from the results of the 
splitting theorem thst significant simplification results from the use 
of three "auxiliary" variablea 4^ , 4^ » *nd A , since each component 
is, in the end, associated with one variable only. Finally, in previous 
works, the fluid is assumed to be a perfect gas with a Prandtl number 
equal to 3/4 , whereas, no such restrictions are imposed here. We now 
briefly examine a few cases where simpler versions of the equations for 
* , * „ , .nd A msy be obtained. 


2.4. Some Particular Cases . 

2.4.1. Prandtl Number ■ Wq/Wq/Pq) : When the bulk viscosity is 

zero, this coincides with situstlons where the Prandtl number is equal 
to 3/4 . Air, which has a Prandtl number equal to .72 , closely ap- 
prosches this case, as long ss its bulk viscosity may be neglected . The 
two diffusion operators on the right-hand side of (2.38) and the left- 
hand side of (2.39) are then identical so that equation (2.39) may be 
used to eliminate 4 in Eq. (2.38). The governing equations reduce 
to: 



1 , Vo .2 

2 >2 + 2 V St 

fl o 3t p o a o 


(2.44) 


at p„ th 



-l)v\ 


(2.45) 


at 


2 -* 

V A 


(2.46) 


The acoustic potential satisfies the modified wave equation (2.44) with 
an "effective" dilatational viscosity Vq^q/Pq . The decoupling is in- 
complete, however, since (2.45) still relates 4 and 4 . As men- 

B tn 

tioned above, this is the situation considered by Wu (1956) and Chu and 
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Kovasznay (1958). 


2,4,2, Non -Ha at -Conducting Viscous Fluid : The Prandtl number of 

Che fluid Is Chen equal co infinity. The encropy is conscanc and che 
Chennai componenc is identically zero. The governing equations, there- 
fore reduce to: 


□ g • - 0 (2.47) 

^ - V Q 7 2 A - 0 (2.48) 


This case was extensively studied by Lagerstrom et al. (1949) and equa- 
tions identical to (2.47) and (2.48) were derived by several authors, in 
particular, Cohen and Tu (1962), Gerlach and Parker (1967), and Scarton 
and Rouleau (1973). 

It is interesting to note that wave phenomena in elastic solids are 
governed by equations which are analogous to (2.47) and (2.48). The dis- 
placement vector 

u * grad $ + curl A (2.49) 


is then shown to satisfy the displacement equations of motion provided 
that the Lame potentials, $(r,t) and A(r,t), are solutions of two 
wave equations with different characteristic velocities, c^ and c T : 

n 1 ^^4 

V •— - o (2.50) 

c, at 


,2 J. _!_»!& 
2-2 
c T at 


(2.51) 


The absence of any diffusion operator or attenuation term in the above 
system Is a consequence of the elasticity of the medium which does not 
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allow for dissipative phenomena. A detailed account of elastodynamic 
theory is given by Achenbach (1973). 

2.4.3. Viscous and Perfectly Heat-Conducting Fluid : The Prandtl 

number of the fluid is assumed to be zero so that there are no tempera- 
ture fluctuations. In such a situation, any wave motion is isothermal. 
Equations (2.38) -(2.40) become: 

□ l*-0 (2.52) 

*th " <V 1)$ a (2 ' 53) 

(2.54) 


2 

where 0 T is the isothermal wave operator defined in (2.43). As seen 
from (2.52) and (2.53), plane longitudinal waves, which, in an inviscid, 
non-heat-conducting medium propagate isentropically with Laplace* s velo- 
city , propagate isothermally in a perfect conductor, with Newton's 
velocity, a^ 'W- 

2.4.4. Barotropic Fluid : A fluid is said to be barotropic when 

k 

the ratio of its specific heats y Q is equal to unity. As mentioned 
by Truesdell (1953), many liquids such as pure water may be assumed to 
be barotropic with a reasonably good approximation. The linearized equa- 


tions of state (2.10)' and (2.11)' now become 
2 . 


P' = a r 


T' = t- 2 - S' 
P 0 


(2.55) 

(2.56) 


More commonly, a fluid is said to be barotropic when its density is a 
function only of pressure. From the expression for Op/5S)p given in 
equation (2.10)', one can immediately establish that these two defini- 
tions are strictly equivalent. 


j 
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and there is no need to define a thermal potential $ h as before since 
the equation of mass and the equation of motion will not contain S' . 
The thermal potential is therefore identically set equal to zero in 
(2.32-2.40), and the entropy fluctuation S' is chosen to characterize 
the thermal component. The governing equations then reduce to: 

V' * grad $ + curl A (2.57) 


^ sa 2 ^ s 
at o at 

T, 


•Vo V *a 


T' 


2 s- 

C 

Pn 


div A = 0 


□sV 1 


-0 

at p 


H - „ 0 V 2 A.° 


(2.58) 

(2.59) 

(2.60) 
(2.61) 

(2.62) 

(2.63) 


Equation (2.62) was obtained by substituting (2.56) into the equation of 
entropy (2.9). 

As discussed previously, the acoustic component is associated with 
pressure, density and longitudinal velocity fluctuations. But in a con- 
trast with the general case, it does not give rise to temperature fluctu- 
ations. The acoustic potential $ is now governed by the modified wave 

Q 

equation (2*61), as in paragraph 2.4.2. 

The thermal component is characterized by entropy and temperature 
perturbations only, and does not give rise to velocity fluctuations. It 
is governed by the classical diffusion equation (2.62). Moreover, 
as shown in (2.61) and (2.62), acoustic and thermal components are com- 
pletely uncoupled, so that an acoustic wave will not generate entropy 
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fluctuations and conversely. This feature persists even in the presence 
of solid surfaces, since it is preserved in the boundary conditions on 
velocity and temperature as seen from (2.57) and (2.59). The analytical 
relationships between thermal and acoustic components are therefore con- 
siderably simplified. 

2.4.5. Inviscid, Non -Heat -Conducting Fluid : In this case it may 

immediately be verified that we recover the basic results mentioned in 
the introduction to this chapter. The acoustic potential is now identi- 
cal to the usual velocity potential $ , and the thermal and viscous 
components are identically zero. Moreover, the state of the fluid is 
governed by the classical wave equation (2.1). 

In the next chapter, we proceed to apply the results of Section 2.3 
to the propagation of small perturbations in cylindrical tubes. 
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III. SOUND PROPAGATION IN A 
TWO-DIMENSIONAL DUCT 

3.1 Introduction . 

The results of the preceding chapter provide a theoretical frame- 
work for the study of the small -amplitude motions of a fluid in a cy- 
lindrical duct. In this part of the investigation, we first derive the 
basic equations governing sound propagation through a duct of constant 
cross-section and then investigate in detail the particular case of a 
two-dimensional duct. 


3.2 Formulation: Parameters of the Problem . 

Consider a harmonic wave of circular frequency u) propagating in 
an infinite tube of constant cross-section S and typical cross-sec- 
tional diameter d (see Figure 1). The fluid is assumed to obey re- 
strictions (a) -(f) discussed in Section 2.2, and in addition, the 
walls of the duct are assumed to be rigid, impermeable, and isothermal. 

The independent variables are non-dimens iona lized as follows: 



(3.1) 


and three non-dimensional parameters are introduced: 



the reduced frequency 


(3.2) 


R 

V 





the “Reynolds number" based on the (3.3) 
speed of sound and the shear vis- 
cosity . 

the "Reynolds number" based on the (3.4) 
speed of sound and the dilatational 
viscosity Tj . 
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The problem therefore depends on the values of five parameters: Vq ; 

P r ; ; R ^ ; and k , which are functions of the state of the medium, 

the dimensions of the tube, and the characteristic frequency of the 
fluctuations. The dependent variables are kept in dimensional form . 
The basic equations of Section 2.3 then take the following form: 

V' * ^ [ grad + grad 4> th + curl A] (3.5) 


with 


v 2 $ 

fit kd a 


Ss 1 


St 



ST 1 = P r R y 
St d 



0 


3* 


th 


(v „- 1)o p 0 st 


divA = 0 


(3.6) 

(3.7) 

(3.8) 

(3.9) 


(3.10) 


and where 4^ , 4) , and A obey the following system: 


[— I 2 4> = k 2 — [ — — v 2 4> ] 

U S a k St L St kR_ V ®th J 


St L St " kR 


S4>. 


th 


St kP R 
r v 


1 v 2 4> 


V 

v o _1 


- - $ 

th kP R a 

r V 


— — A = 0 

3t kR A U 
V 


at 


□s { * or *th 5 ■ kpV ' 7 bit* s ° r *th 1 


r v 


(3.11) 

( 3 . 12 ) 

(3.13) 

(3.14) 
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where 



k 


2 a 2 k 

at 2 \ 






(3.15) 

(3.16) 


In the above relations, the primes on x' . y 1 , z 1 . t 1 , have 
been omitted for convenience . The system (3. 10) -(3. 14) is to be sup- 
plemented with the following boundary conditions at the walls: 



on duct boundary walls. 


(3.17) 


We wish to determine the potentials $ a (r,t) , $ t ^(r,t) , and A(r,t) 
which satisfy (3. 10) -(3. 14) and the associated boundary conditions (3.17). 
The previous considerations define an eigenvalue problem which is 

— - —4 -4 -4 

more easily handled if 9 (r,t) , 9 (r,t) , and A(r,t) are assumed 

a tn 

to be travelling waves of the form 

i(t-fl z) 

♦ (r,t) = cp (x,y)e (3.18) 

a a 




i(t 


cp th (x,y)e 


• z) 
z 


(3.19) 


_ _ i(t-fi z) 

A(r,t) =G (x,y)e 


(3.20) 


Note that in these relations, both z and t are non-dimensionalized. 
The circular frequency o> is embedded in the non-dimensional time t , 
and the complex wave number is non-dimensionalized with respect to 

the duct width. Equations (3.18)-(3.20) physically correspond to stand- 
ing waves in the cross-section of the tube and travelling waves along 

the duct axis z. . The real part of 8 denotes the actual propagation 

z 
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wave number, while the imaginary part represents the attenuation rate. 
The complex quantity j3 , and the functions cp (x,y) , cp . (x,y) , and 

2 Q tn 

G(x,y) are to be determined as a function of the five previously de- 
fined parameters. We now proceed to derive the governing equations for 
cp fl , cp th , and 2 . 

Substitution of (3.18) or (3.19) into the fourth order partial dif- 
ferential equation (3.14) yields the following relation: 


iv 0 k ik 


7 1 ^a or ^ " C ikP r R i/ 1 + iv.k/R^ 


where 


o o , 1+iv k/P R +ik/R , 

+ 2« 2 ]7 1 t» a or<p tli ) + CS 2 +lUP r R v 71 S, 2 (3.21) 


ik 3 P R 


0 7J 


- 5 

C X hx 



(3.22) 


Solutions of (3.21) may be obtained in the form of a combination of solu- 
tions of the following Helmholtz equations: 


2 2 
V®ao° r Vo 1 + “<>.<, ° r<p tJ 

= 0 



(3.23) 

V l ( ‘ p s l or<l> thl ,+0! l t ' p 8 l or<p thl 1 

= 0 



(3.24) 

2 2 

where and Oc ^ are the roots of the 

algebraic equation in 

a 2 , 

2 

obtained from (3.21) by changing and 


2 

into -a 

and 

a 4 > 

respectively. The acoustic potential cp 

a 

will 

therefore be 

the 

sum of 

cp and cp , , solutions of (3.23) and (3. 

ao a l 

24). 

Similarly, 

’th 

will 
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be the sum of cp . and cp . , . 

T tho T thl 

It Is a straight-forward procedure to derive expressions for Qt ( 
2 

and from (3.21). The results are as follows: 


a o 2 • »o ■ ^ < 3 - 25 > > “i 2 ■ ^-K 2 


(3.26) 


with 


, -ikP R V_k . 

*o,l " 2 ( l+iV^k/R^) [1+1 PT + 1 ^ 



1 - 


4ik 


i-nypk/R^ 


P r \ ( l+iV Q k/ P r R y + ik/ O 2 J 


(3.27) 


where the o and 1 subscripts correspond to the - and + signs 
respectively. 

In addition, the acoustic and thermal potentials cp^ and cp t ^ are 
coupled through (3.11) and (3.12). When (3.18) and (3.19) are substi- 
tuted into (3.12) and use is made of (3.23) and (3.24), one easily ar- 
rives at the following relations: 


P n + ikP R 

Cf ao (x,y ) = *~2 W X ’ y) 

<V 1)/3 o 2 


(3.28) 


cp a i (x , y) 


0. +ikP R 

- 1 ¥ W x ’ y) 

(y 0 -l)^ 2 thl 


(3.29) 


One then may check that the second coupling relation, (3.11), is also 
satisfied . 

In a similar fashion, when use is made of (3.20) in (3.13) and 
(3.10), one obtains the governing equation for the vector potential: 
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(3.30) 

(3.31) 

(3.32) 

with the following restriction on the possible solutions 

V 5 i- ifi z a * < 3 - 33 > 

In the above coupling relation, the vector potential 2 has been de- 
composed into an axial component 2^ and a component in the cross -sec- 
tional plane ♦ A similar notation will subsequently be used for 
the velocity vector V' . 

Finally, the boundary conditions (3.17) may be written in terms of 
the auxiliary variables as follows: 


where 


and 


VjS+a 2 2 2 - 0 


„ 2 * 2 * 2 

a 2 ^2 -K 


ft. * -ikR 
2 v 


r dv{ v x,y) = v (cp + cp ) + V. xG - i/3 e X G = 0 
I 1 a tn 1 z zzl 


dV 


' (x , y) = -i/J (cp + <p . )e + V x2 = 0 ° n ^j Ct 
z ,J/ z a th z 1 1 boundary 


walls. 


V. 



T'(x.y) = cp th 


0 


(3.24) 

(3.35) 


(3.36) 


Discussion : For a cylindrical duct of arbitrary cross-section and given 

values of the parameters Vq , , R^ , and k , the problem has 

been reduced to the determination of a complex wave number )3 and three 
unknown functions cp fl , <p h , and 2 . The general analytical form of 
the auxiliary functions is given by equations (3.23), (3.24), and (3.30) 
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which have to be supplemented with the coupling relations (3.28) , (3.29) , 

and (3.33). Enforcement of the boundary conditions (3.34) - (3.36) then 

determines the remaining unknown /9 or arbitrary cross-sectional 

2 

shapes, it is unfortunately impossible to find separable solutions of 
equations (3.23), (3.24), and (3.30), compatible with the boundary con- 
ditions at the walls, and leading to a dispersion relation for /9 ex- 

z 

pressible in terms of a finite number of transcendental functions. It 
was shown in Section 2.4.2 that the analytical formulation of wave pro- 
pagation problems in elastic solids bore a close resemblance with the 
model presented in Chapter II. The same mathematical difficulties have 
indeed been encountered by researchers interested in the propagation of 
elastic fluctuations in solid wave guides. A review of this subject is 
given by Meeker and Meitzler (1964) and Achenbach (1973). However, in 
the case of a two-dimensional or a circular tube, one may arrive at an 
exact dispersion relation in terras of a finite number of tabulated 
functions . 

We have chosen to treat here the propagation of waves between two 
parallel infinite plates, i.e., through a two-dimensional duct. Many 
similarities will be shown to exist between the families of modes in a 
two-dimensional duct and a circular tube, and identical methods may be 
used in both instances. Since we are mainly interested in bringing out 
the main physical features of the problem, and wish to minimize analyti- 
cal complications, we devote the major part of this study to a detailed 
investigation of wave propagation in a two-dimensional duct. 

3.3 Two-Dimensional Duct : 

Let us consider a two-dimensional duct of width d as shown on 
Figure 2 and a harmonic wave of frequency u> propagating in the positive 
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z direction, under the seme assumptions as in the previous section. 

In analogy with equations (3.18) - (3.20), we further restrict the 
functional forms of $ , ♦ ^ , and A to the following: 


* (r,t) - « (x)« z 

a fl 


(3.37) 


*th (t>t) * 'W 0 * 


i ( t z ) 


(3.38) 


A(r,t) - G(x)e 


i(t-fiz) 


(3.39) 


Note that cp , cp . , and Ct are assumed to be independent of the space 

a tn 

coordinate y . In other words , the wave motion Is the same In all planes 
normal to the y-axis. However, since we allow the vector potential to 
have x- and z -components, the motion is not restricted to the y-plane, 
and the velocity field may admit a non-zero component along the y-axls . 

The governing equations (3*23), (3.24), and (3.30) simply become: 


— - {cp or cp . } +0L. {cp or cp . } « 0 

, 2 ao tho 0 LY ao ^tho 
dx 


( 3 . 40 ) 


“2 Kl or <P thl ,+ “l t<p .l or "’thl 1 ■ 0 
dx 


( 3 . 41 ) 


4 +a 2 2 * o 

dx 1 


( 3 . 42 ) 


so that the general form of <P. h («) is: 


co . (x) «A cos a x + Bsina x+C cos a, x + D sina,x 
th o o 1 1 


( 3 . 43 ) 


and the corresponding form of cp (x) derived from the coupling relations 

a 

(3.28) and (3.29) is: 
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9.(x) 




' .>Vj 






r 

► 


& +ikP R 

[A cos a.x + B sinO A x] 

<v„-»A,‘ 0 0 


0.-HkP_R 


+ -- 


«V l >* 


[ C cos a. x + D s In a. x] 

A A 1 1 


(3.44) 


Similarly, use of (3.42) and the condition (3.33) yields the following 
expressions for the components of the viscous potential: 


I 


G (x) • Ecoia x + Fjlna x 

X 4 4m 

(3.45) 

G^(x) « G cosa^x + H sinC^x 

(3.46) 

^2 


G (x) * (- F cosO^x +E sinO^x) 

(3.47) 


z 


Finally, the boundary conditions (3.34) - (3.36) become: 


v; ( x> « s «p ii «r <p th ) +i$a y - o 


dG 

z 

dx 


v;.(x) « - -rr - i/8 C - o 


Z X 


dG 


at x ■ + \ 


V'(x) oc -if} (cp +Cp ) + -r* ■ 0 
z z a tn dx 

T * (x) « 9 th - 0 


(3.48) 

(3.49) 

(3.50) 

(3.51) 


and substitution of (3.43) - (3.47) into (3.48) - (3.51) yields the fol- 
lowing set of eight homogeneous linear equations for the eight unknown 
constants: A, B, C, D, E, F, G, and H: 


W +lkP - R 


r V 


[ -A sin “ + B cos 1 +a. 


y^3, 2 +ikP_R 


+/\ n L -** »■*■** * * tj vwo 0 j 

° <V‘>0o 


a a a a 

[ -C sin-j-+D cos — ] + ifl [G cos — + H sin — ] * 0 


(3.52) 
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yJ3 n 2 +ikPji t ,_ % a n ^ Vl +ikP A 

a n — 2 ^ [A sin — +B cos — ] +aj 


‘0 






Of. a i tt 2 a 2 i 

[C sin-r* + D cos -r* ] + ip [G cos-r- - H sin-r* ] - 0 


(3.53) 


yfl 2 +ikP R QC Ot V/A +ikP R 

•1/S YCAcosf • 

2 (V.-Dft, 2 2 


a, a , «, “ 9l 

[C cos-^+D sin-^] +a 9 [-G sin-j-+Hcos — j = 0 (3.54) 


v^„ 2 +ikPfi, r « n a o, V i l +ikP r K l/ . 

-i/3 _0_0 E-K[Acos^-Bsin T ] - TT 

2 <V 1)(J o 2 2 <V 1)(l i 


a. a, a 2 a 2 


[C cos sin-j 2 3 +0C 9 [G sin -y- -Hi cos — 1 * 0 (3.55) 


a a Q p 1 Qfj 

cos +B sin — + C cos — + D sin-y = 0 


(3.56) 


a & Q a l 

A cos -v-Bsin-^- + C cos — - D sin— = 0 


(3.57) 


a o' 

2 2 

E cos — + F sin— = 0 


(3.58) 


a 2 “2 
E cos — -F sin — = 0 


(3.59) 


By elementary manipulations such as addition or subtraction of consecu- 
tive equations the above system may be immediately expanded into the four 
following subsystems. 

System 1 : 


“2 

E cos — = 0 


(3.60) 
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System II : 


“2 

F sin — = 0 


(3.61) 


System III : 


v /3 2 +ikP R a y 8, 2 +ikP R a, a 0 

a r^Asin-^+a — sin-r^- ij8 Hsin-^=0 

<V l) V (y 0 -DA 2 2 2 2 

(3.62) 


.* ^ V... 

2 (Vn-DA 


r 0 /r X) 

oc 


o « ~ o C cos — +G£ 0 H cos “*t~ — 0 

2 %-Dfii 2 2 2 


a. 


A cos + C cos ~y - 0 


(3.63) 

(3.64) 


System IV : 


0L r 


yj3 2 +ikP R 0c n v 8 2 +ikP R a, 

a V^Bcos-^+a -2-i ^Dcos-^+iiS G cos 0 

<V 1)/3 o (v n V 2 2 2 

(3.65) 


•ifi 


V» wkp A r . “o . “i “ 2 

r— Bsm — - x^ r^Dsin — - a.Gsin — = 0 


2 O'o- 1 >«o 


2 <r o- 1)3 i 


2 2 


a o a i 

B sin-j- + D sin— = 0 


(3.66) 


(3.67) 


In order for these subsystems to admit non-trivial solutions, the deter- 
minants formed with the coefficients of their respective variables must 
be set equal to zero. Those conditions yield four compatibility condi- 
tions or so-called dispersion relations which may be written as follows: 
Dispersion Relation I : 

Oi, 


2 „ 
cos — = 0 


(3.68) 
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Dispersion Relation II : 



(3.69) 


Dispersion Relation III : 


2 1 i - 2 “o ^1 °^2 k 2 ^o k 

k j co, y ,ln 7*<“T ■ ff* °o a ! 

8 o r v 


a a a 2ivk a a a 

.0 1 2 . k ~VV „ “o . 1 2 _ 

sin y cos y cos y-^- — ^ cosysinycosy = 0 

*1 r V (3.70) 


Dispersion Relation IV : 


2 1 l 2 % ^1 ^2 

k( 7T‘T2 ) ^ T sin -cosy +(- 5 -—)^ 

*0 $ 0 r V 
“O “l “2 k 2 iV 0 k “o “l “2 

cosysinysiny-(- ? -— ) ttl a 2 sin ^cosysiny = 0 

p 1 xv 

(3.71) 


Discussion : As seen from equations (3.60) - (3.67), each subsystem 

pertains to distinct groups of unknown constants. Consequently, if any 
one of the dispersion relations (I) to (IV) is satisfied, there exists 
a non-zero solution for the unknowns of the corresponding subsystem. 

At the same time, since the remaining determinants then take non-zero 
values, the only possible solution for the other subsystems is zero. 
Hence, each set (I) to (IV) and its dispersion relation may be inves- 
tigated separately, the constants associated with the other sets being 
identically zero. Even though a detailed description of the four cor- 
responding types of wave motion will be given in Section 3.4 and Chapter 
IV, we take advantage of this mathematical feature to briefly outline 
their main characteristics. 
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Let us first indroduce the following definitions: A wave motion 

will be called symmetric (respectively, antisymmetric) if its velocity 
is symmetric (respectively, antisymmetric) with respect to the median 
y-z plane. It will be referred to as inplane (respectively, antiplane) 
when its velocity vector lies in (respectively, is perpendicular to) the 
x-z plane. From the above discussion, and the form of equations (3.60)- 
(3.61), we infer that, in sets (I) and (II), all constants are zero with 
the exception of E or F . As seen from (3.43) - (3.47), the acoustic 
and thermal components are then identically zero, and the vector poten- 
tial lies in the x-z plane. The velocity vector is therefore purely 
rotational and its only non-zero component is V given by (3.49). 

y 

Furthermore, F is zero in System (I) so that is then an even 

function of x . Likewise, E is zero in System (II), so that is 

then an odd function of x . We may conclude that sets (I) and (II) de- 
scribe the behavior of antiplane symmetric and antisymmetric velocity 
fluctuations. In a similar fashion, when one considers set (III) (re- 
spectively, (IV)), B, D, E, F, G, (respectively A, C, E, F, H) are zero. 
Upon examination of (3.43) - (3.50), it is clear that both sets describe 
inplane wave motions, and that (III) and (IV) are associated with sym- 
metric and antisymmetric velocity fields respectively. However, in 
contrast with the antiplane motions, acoustic, thermal, and viscous 
components are now inherently coupled and lead to complicated wave con- 
figurations. The initial problem of wave propagation between two infi- 
nite parallel walls has thus been reduced to four simpler subproblems, 
namely, inplane or antiplane, symmetric or antisymmetric wave motions, 
each one being governed by one of the Systems (I) to (IV) and its cor- 
responding dispersion relation. 
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In equations (3.68) - (3.71), the complex quantities a Q , and 
U may be expressed in terms of the complex wave number fi through 
the use of (3.25), (3.26), and (3.31). Each one of the dispersion re- 
lations is then an equation for /3 , and its solutions are to be deter- 

z 

mined as a function of the given non-dimensional parameters , P r , 

, R^ and k . Each solution defines a mode of propagation, char- 
acterized by specific variations for cp (x) , cp , (x) and G(x) . Such 

8 til 

characteristic variations constitute a mode shape, and they may be ob- 
tained exactly by solving the appropriate subsystem (I) - (IV) for the 
corresponding unknown constants. Equations (3.68) - (3.69) are trivial 
and the next section will examine their solutions and associated mode 
shapes. On the other hand, the dispersion relations (3.70) - (3.71) 
pertaining to the inplane modes, are transcendental equations for , 
and cannot be solved exactly. In Chapter IV, we will propose a pertur- 
bation scheme, whereby approximate solutions may be obtained. 

We have previously pointed out the analogy between wave propagation 
in viscous fluids and in elastic solids. It is indeed striking to note 
that wave motion in an elastic layer is governed by equations which are 
very similar to the subsystems developed here. Meeker and Meitzler 
(1964), in particular, assumed the Lame potentials to have a form anal- 
ogous to the expressions chosen in (3.37) - (3.39), and were thus led 
to classify the possible solutions in terms of four families of waves, 
namely, symmetric and antisymmetric SH (Shear-Horizontal) waves, and 
longitudinal and flexural plane strain waves. The dispersion relations 
describing the latter two families, or so-called Rayleigh-Lamb equations 
are found to present essentially the same basic features as relations 
(3.70) and (3.71) when heat flow is assumed to be zero. In such a case, 
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the last term in relations (III) and (IV) disappears. In the next sec- 
tion, we discuss the characteristics of the symmetric and antisymmetric 
antiplane wave motions. 

3.4 Symmetric and Antisymmetric Antiplane Wave Motion . 

The solutions of dispersion relation (I) are immediately given by: 

a = (2n+l)TT (3.72) 

when n is a zero or positive integer. Hence, there exists an infinity 

of symmetric antiplane SA-modes. Their shape in any cross-sectional 

plane is characterized by a specific value of the coefficient OL ^ . The 

propagation wave number attached to each mode is obtained by substituting 

(3.72) into (3.31). In solving the resulting equation for )3 , we se- 

z 

lect the complex root which pertains to waves attenuating in the posi- 
tive z-direction, i.e. of negative imaginary part, the other root being 

* 

associated with waves attenuating in the negative z-direction. Such a 
choice does not restrict the scope of the study, since the propagation 
characteristics are obviously independent of the direction along the 
z-axis. The wave number is then given by the following relation: 



S~2 

-i^ 


7 / 


( ( 2n+l )TT) 4 +(kR^) 2 - ( ( 2n+l)n) 2 


( ( 2n+l ) tt) 4 +( kR y ) 2 +( ( 2n+l ) tt) 2 


(3.73) 


We have avoided here using a radiation condition based on the sign of 
the phase velocity, i.e. the direction of propagation of the wavefronts, 
since, in Chapter IV, we will encounter backward -propagating waves which 
decay in the positive z-direction. 
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As discussed in the preceding section, E is the only non-zero constent, 
so that from (3.43) - (3.47), acoustic and thermal potentials are iden- 
tically zero, and the viscous potential lies in the x-z plane with the 
following components: 


SA 

G (x) ■ E cos ((2n+l)Trx) 
xn 


(3.74) 


/-SA. . 

G (x) - 0 

yn 


,-SA, s (2n+l)rr i „ , ,, n ,, x v 

G zn ( x ) * — E sin ( (2n+l)nx) 

^n 

the rotational velocity field is then such that: 

-» SA -a 

V* (x) * Jr E cos ((2n+l)nx)e 

n di3 SA y 

zn 


(3.75) 


(3.76) 


(3.77) 


All other physical variables are identically zero. 

The velocity profile of each mode will be normalized by the follow- 
ing condition 


V^ n (0) - 1 (3.78) 

which determines the value of the arbitrary constant E . The final form 
of the mode shape is: 

SA 

V yn ^ * cos(2n+l)nx (3.79) 

The characteristics of the symmetric antiplane mode have therefore been 
derived in a straight-forward manner due to the extreme simplicity of 
the associated dispersion relation. They are purely transversal rota- 
tional waves, the velocity fluctuations being perpendicular to the plane 
of propagation. An interesting interpretation of equation (3.79) follows 
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from the decomposition of the cosine function into complex exponentials. 
The velocity may then be written as follows: 


y |SA 

yn 


(x,z,t) 



i[t + (2n+l)rrx - 
e 


Re z] 
zn 


- Im 8 SA z i[t - (2n+l)nx - Re z] 

+ £e 2n -e zn (3*80) 

SA 

Each mode may thus be considered as the superposition of two plane 

transversal waves travelling in symmetric directions with respect to the 
median yz -plane and, which* upon multiple reflections at the boundaries, 
give rise to a standing wave pattern in the cross-sectional plane, and 
a propagating wave along the z-axis. 

Analogous results may be obtained for the antisymmetric antiplane 
AA -modes. Their complex wave number is then given by: 



1 




(2nTT) 4 +(kR y ) 2 -(2nrr) 2 


(2nTT) 4 +(kR y ) 2 +(2rm) ‘ 


(3.81) 


where n is a positive integer, and the corresponding mode shape is: 


V ,AA (x) = sin2nTTx 
yn 


(3.82) 


The magnitude of the imaginary part of the complex wave number de- 
fines the non-dimensional attenuation rate per diameter along the duct 


* 

The integer n cannot be zero, since such a value leads to a mode 
shape which is identically zero. 
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ATTENUATION RATE 



Figure 3. Symmetric Antiplane Modes. Attenuation rate 
versus reduced frequency for R v » 2.35 x 10^, 
d * 10"^, medium is air at 15°C, 1 atm. 
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PHASE VELOCITY 



Figure 4* Symmetric Antiplane Modes. Phase velocity 

versus reduced frequency* Same values of the 
parameters as in Figure 3. 
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axis, whereas the real part Is the actual propagation wave number. The 
phase velocity non-dlmensionalized with respect to the speed of sound 
may be obtained by the formula: 




(3.83) 


On Figures 3 and 4, the attenuation rate and the phase velocity of the 
first three symmetric antiplane modes are represented as a function of 
the reduced frequency k , for a given set of values of the other four 
parameters. These curves describe the dispersive properties of the 
waves, and will be used extensively in the next chapter. As seen from 
these plots, one may distinguish two ranges of reduced frequencies: 
the low-frequency range where both attenuation rate and phase velocity 
are constant and independent of k , and the high-frequency range where 
they increase as the square root of k . These two regions are sepa- 
rated by a transition region. Further details on this question will be 
given in Chapter IV. At any rate, it is apparent from Figure 3, that 
antiplane modes are highly attenuated at all frequencies. 

The antisymmetric antiplane modes have strictly similar character- 
istics, and since this study focuses on the symmetric modes, we do not 
need to elaborate on their properties. One may finally note that equa- 
tions (3.73) and (3.81) may be merged into a single relation in terms 
of a new index N , equal to 2n + 1 for the symmetric modes and 2n 
for the antisymmetric modes. The eigenvalues of the antiplane modes 
are then given as follows: 

? ! 
zn 

(3.84) 




+(kR y ) -(NTT)' 


- i 


+ < k V‘ 


+ (Ntt) ‘ 
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where 


N ■ 1, 3, 5, ... for the SA inodes 
and 

N » 2, 4, 6, ... for the AA modes 
3.5 Conclusion . 

For the purpose of Illustrating the basic method of solution, we 
have treated in detail the simple case of the antiplane modes. In this 
instance, exact solutions of the dispersion relation were immediately 
obtained throughout the entire reduced frequency domain. Unfortunately, 
dispersion relations (III) and (IV) are much more complicated, and re- 
quire the use of some approximation scheme, if any analytical results 
are to be obtained. Chapter IV is, therefore, devoted to a perturbation 
study of the inplane wave motions, with particular emphasis on the sym- 
metric inplane modes. 
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IV. PERTURBATION STUDY OF THE INPLANE MODES 


4.1 Introduction . 

In was shown In Section 3.3 that the eigenvalues associated with the 
symmetric (respectively, antisymmetric) inplane modes are solutions of 
dispersion relation (III) (respectively, (IV)), and that the corresponding 
eigenfunctions are given by (3.43) - (3.47) where the coefficients A , 

C , H (respectively, B , D , G ) satisfy System (III) (respectively, (IV)), 
and E and F are zero. This part of the investigation will be mainly 
concerned with the derivation of approximate solutions of the transcendental 
dispersion relations (111) and (IV). From these solutions, a detailed pic- 
ture of the possible inplane wave motions will emerge. 

Before undertaking such a study, we note that, when dispersion rela- 
tion (III) is satisfied for a given eigenvalue /3 , System (III) admits 

z 

an infinity of non-zero solutions which may be obtained by solving equations 
(3.62) and (3.64) for C and H . We have then the following relations: 


C 


-os a Q /2 
cos 0^/2 


(4.1) 


II 


kP rV 1 

V 1 V 2 


_J_.fr co> V 2 , 

a 2 } a cos a /2 


(4.2) 


where A is still an arbitrary constant and will be determined later by 
imposing a suitable normalization condition. Substitution of the above 


Tlw. .notion is therefore purely two-dimensional and the vector potential 
reduces to a stream function A . In an earlier report, Huerre and 
Karamcheti (1973) a priori considered such a stream function, and thus 
restricted their investigation to inplane motions. The slightly more 
general approach used in this study presents the advantage of also re- 
vealing the existence of antiplane wave motions. 
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values for C and H Into (3.43), (3.44), and (3.46) leads to the expres- 
sions for cp (x) , q> . (x) and C (x) . The corresponding variations of 
a tn y 

tha phyalcal variables are then given by (3.5) - (3.9). The resulting 
symmetric inplsne mode shapes are shown in Table 1. We must emphasize 
that such mode shapea are functions of the specific wave number A and 
that they cannot be determined completely until the dispersion relation 
has been solved. However, when the eigenvalue is known, the associated 
mode shape will be obtained immediately by substitution into the relations 
of Table I. The same procedure may be used to determine the antisymmetric 
mode shapes given in Taule II. We now proceed to an in-depth analysis of 
the inplane eigenvalues. 


4.2 Methodology - Preliminary Assumptions . 

In order to display the common features of the dispersion relations 
pertaining to the Inplane wave motions, wc divide (3.70) by o^cosa^ . 
cosOfj/2 cosOt 2 /2 and (3.71) by a 2 sina Q /2 sinO^/2 sina^/2 , and recast 
them into the following single formula: 


2 2 Q a 2tvk a 

, k V . Zl f t*n i _2. ( ±_ 7 0 r tan A 

“ R 2 } 0C, l cotan J 2 K R 2 “o l cotan J 2 

P 0 


0 , 


0r 


r v 


, k 2 iy 0 k . t tan A 

' 2 ' P R )fl i cotan 2 

P, TV 


(4 . 3a ,b) 


where the upper and lower trigonometric functions refer to the symmetric 
(4.3a) and antisynmetric (4.3b) dispersion relations, respectively. In 
moat instances, we will only need to consider the above expressions. How- 
ever, the normalization which has just been performed may have restricted 
the number of possible solutions of the initial relations (III) and (IV), 
and in a few cases, we will find it more convenient to study the alternate 
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TABI£ II. Antisymmetric Inplane Mode Shapes. The 

constant B has been replaced by 
iP R a 

B * = (yf-m B sin T • 


o" 
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set of relations: 


o 

oo B a a 2ivk 

, k Jtv P z rcotan^ Q rcotan 1 . , k 0 ^ . 

{ Q 2 * 2 ; a o 1 tan j 2 1 tan J 2 V"PR j 0 

B o 8 i 2 p o r v 

ft ft 2 iV k 0L OL 

r cotani _1 rcotan-i _2 , k 0 n .. fcotani _0 rcotan-i _2 _ Q 

1 tan J 2 1 tan J 2 “ ( B 2 'P R / “l 1 tan J 2 1 tan j 2 

Pi r j/ 

1 (4.3c,d) 

with the same conventions for the trigonometric functions as in (4.3a,b). 

This latter form of the dispersion relations results from the divisions 

«Q a i a 2 a 0 a i a 2 

of (3.70) by a 2 sin— sin — sin— and (3.71) by a.^ cos — cos -y cos — . 

Together with the first normalized set, it will ensure that all possible 

solutions of the original equations are examined. 

Expression (4.3a) is strikingly similar to the dispersion relation 
first derived by Kirchhoff (1868), and describing the axisymmetric modes 
of a duct of circular cross-section. To obtain the latter, one simply 
replaces the trigonometric tangent in (4.3a) by the ratio of Bessel func- 
tions, J^/ Jq . We, therefore, expect the symmetric inplane modes to be 
closely related to the axisymmetric modes of the circular geometry. 

As mentioned in Chapter III, both equations are transcendental and 
cannot be solved exactly. This constitutes the major obstacle of the 
present investigation. The problem may be approached from two basically 
different ways. A numerical scheme can be developed to isolate the solu- 
tions of the equations. Such a method was implemented by Shields et al. 
(1965) and Tijdeman (1975) to determine the characteristics of the fun- 
damental zeroth-pressure mode, and by Scarton and Rouleau (1973) to 
study all the axisymmetric eigenfunctions in the case of zero heat-con- 
duction. The analysis is then considerably simplified by the presence 
of only two terms in the dispersion relations. In a second approach sug- 
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gested by Kirchhoff (1868) and Rayleigh (1877), and adopted by a large 
number of investigators to be mentioned later in the course of this 
chapter, one attempts to obtain approximate solutions of the dispersion 
relations analytically . Applications of this type of procedure have 
traditionally been restricted to the determination of the characteristics 
of the pressure modes in the acoustic boundary layer approximation. The 
pressure modes, in the limit of zero boundary layer thickness, reduce to 
the familiar acoustic modes encountered in inviscid propagation problems, 
so that the first terms in the approximations are then readily available. 
In this study, we develop a methodology whereby perturbation expansions 
of the complex wave number may be extracted from the dispersion relation 
for all the eigenfunctions of the problem, in as many ranges of the para- 
meters as possible. The main advantage of this approach over the purely 
numerical one is that it exhibits in a compact mathematical form the es- 
sential physical features of the inplane modes as well as their dependence 
on a few non-dimensional parameters. Before describing the details of the 
method, we draw the consequences of some of the assumptions made in Chap- 
ter II. 

Preliminary Assumptions . Let us seek, in the particular case of a 
gaseous medium, an estimate of the non-dimensional quantities and 

k/Ry , where k and R^ are defined by equations (3.2) and (3.3). From 
elementary kinetic theory considerations, we know that the kinematic vis- 
cosity Vq is such that 


v o - V 


(4.4) 


where l is the mean free path. This relation yields the following ap- 
proximation for and k/R^ : 
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(4.5); 


(4.6) 



_k_ oj£ l 

R “ a. ~ X 
V o 


where X is the wavelength associated with the circular frequency co . 
Hence, and k/R^ are respectively a measure of the ratio of duct 

width to mean free path and mean free path to wavelength. In order to be 
able to write the basic equations of Fluid Mechanics, we assumed in Sec- 
tion 2.2 that the fluid was a continuum. Consequently, in the case of a 
gaseous medium, the width of the duct and the characteristic wavelength 
of the fluctuations X are both assumed to be much larger than the mean 
free path, so that we have: 

R » 1 ; jr«l (4. 7a ,b) 

V K , 


Without any additional assumptions other that those stated in Chapter II, 
we may then consider (4.7a) and (4.7b) to be satisfied. Hence, the pre- 
sent investigation is restricted to values of the frequency parameter k 
such that 


0 < k « R^ (4.8) 

For air at normal pressure and temperature, the above inequalities imply 
the following limits on the duct diameter and frequency of the wave: 

d » 5XlO“ 6 em (4.9) ; f « 10 9 Hz (4.10) 

When we examine the corresponding limits for water under the same 
conditions , 

d » 8X10' 8 cm (4.11); f « 3xl0 U Hz (4.12) 

it is clear that, in all practical situations, 1/R^ ant * k/ R y may be 
considered as very small quantities. The same conclusions may be drawn 
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for 


1/R^ and k/R ^ since the bulk viscosity is never much larger than 
the shear viscosity. Furthermore, and are of the same order of 

magnitude. 

As a result of this discussion, the complicated expression (3.27), 

2 2 

which defines $ Q and ^ can be considerably simplified by expansion 
in powers of k/R y to yield the following: 


o 9 y.-l R , 

r Tf v 


(4.13) 


V ' - ikP rV 1+i<y 0- 1>< f-R li) f :l 

r r\ V 


(4.14) 


When such expressions are substituted into relations (4.3a,d) we obtain: 

P r V a „ P r <p r S,\ ' 

f tan 4/_L ^ jLi/v r tan , 

^ cotan 2 +[ “ i( P r ' R^ R y ■* ^cotan 3 2 + 




and 


P r VV P r t VV “2 ' 

rcotani^o rcotan-> a i r 1 j. k, 

1 tan j 2 1 tan j T + [l " i( F R } T 1 % * 

r T) V 

rcotan-,^1 rcotani ^2 n k r . . , 1 j k, 

tan 2 ** tan ^ 2 + P R [ 1 + i(p - R ) R J 

r v x 7] y 

rcotan-i 0 rcotani 2 
a i { tan } T tan 3 T“ 0 


(4.15a,b) 


(4.15c,d) 


where 
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VVVt 


£ 

“I 2 - -itttR^i+Kyo-u^-^r 3 -^, 


“2 2 • - ik \-^ 2 


(4.16) 

(4.17) 

(4.18) 


and the coefficients of each term in (4.15a,d) have been expanded to 

k ^ 

order (— ) included. 

K V 

Methodology . Even though the previous comments have led to a sub- 
stantial simplification of the dispersion relations, as shown in equations 
(4.15a,d), they do not provide a method of solution. An important indi- 
cation as to how to approach this problem, may be seen in a short exami- 
nation of the physics of wave phenomena in viscous, heat-conducting and 
compressible fluids. As mentioned briefly in Section 1.2, Chu and Ko- 
vasznay (1958) distinguish three main types of fluctuations in such a 
fluid: pressure fluctuations corresponding to the propagation of sound 

waves in an inviscid fluid, vorticity fluctuations related to the diffu- 
sion of vorticity perturbations in a viscous medium, and entropy fluctu- 
ations related to the diffusion of so-called ,r hot spots" in a heat-con- 
ducting medium. In this study, we therefore expect to encounter three 
families of pressure-, entropy-, and vorticity-dominated modes, and we 
are faced with the problem of finding a method which enables us to deter- 
mine their respective characteristics. In inviscid wave propagation, the 
first and last terms in (4.15a,b) are identically zero, and the disper- 
sion relation is immediately solved for , to yield the pressure or 

acoustic modes. In the case of a viscous and heat-conducting medium, 
these terms are no longer negligible. They may be recast solely in terms 
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of 0 £q through the use of (4.16), (4.17), and (4.18), and the resulting 
relation is still to be solved for , to yield the pressure-dominated 
modes. Similarly, the characteristics of the entropy- (respectively, 
vorticity-) dominated modes will be found by recasting the dispersion re- 
lation in terms of 0^ (respectively, ) and solving for 0^ (respec- 
tively, a 2 ). The main advantages of such a procedure, as compared to a 
straightforward solution in terms of the complex wave number /3 , will 
become apparent in the next few sections. 

We are still confronted with the task of expanding the dispersion 
relation, and the unknown; , 0^ , or , in terms of a suitable small 
parameter. Three such parameters may be defined, and they correspond to 
three ranges of frequencies, or equivalently three ranges of duct widths: 
The High-Frequency-Wide -Tube Range . 

Where the reduced frequency k is such that: 

r«k«R„ 1/3 (4.19) 

K .. v 


and the small parameter is defined as: 


€ 



(4.20) 


This is the familiar acoustic boundary layer approximation. 
The Low-Frequency-Narrow-Tube Range . 

Where the reduced frequency k is such that: 


k 



(4.21) 


and the small parameter is defined as 


€ = kR 

V 


(4.22) 
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The Very-High -Frequency-Very -Wide -Tube Range . 

Where the reduced frequency k is such that 


R « k « R 
V V 

and the small parameter is defined as 



(4.23) 


(4.24) 


These ranges were first suggested, in a somewhat different form, by Weston 
(1953a) in an analytical study of the zeroth-pressure mode. A full justi- 
fication of such definitions will arise from the detailed analysis of the 
next three sections. 

The methodology which has just been outlined is now applied to the 
determination of the characteristics of the pressure-, entropy-, and vor- 
ticity -dominated modes in the three ranges of parameters defined above. 


4.3 High Frequency - Wide Tube Range . 

We showed in the preceding section that, in all practical situations, 
1/R is a very small parameter. It is, therefore, legitimate to seek a 
formal expansionof a Q (respectively, 0^ , CC 2 ) in powers of 1/R y . 

The range of validity of such expansions will be examined a posteriori, 
by requiring that the ratio of two successive terms be smaller than unity. 

4.3.1 Pressure -Dominated Modes . We assume the unknown 0 £q to be 
of order (1/R )^ or of higher order, and wish to expand the dispersion 
relation (4.15a,b) to order 1/R^ , inclusively. From (4.16) - (4.18), 
the coefficients j3 , a, , and a may be approximated by the follow- 
ing expressions: 


8 

z 


2 


(k 2 -a r 


-)[ 1 - 


ik 


k 2 -a 


- (-^— + - if ) — ] 
2 ' P R ; R J 
r T) V 


(4.25) 
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n 4 \ / kp r R t/ r 1 4.4 a 0 ~ k , . V 1 /_L ilfv Jli /a 0 fi\ 

(1-1) >y 2 ^ 1+i 2kP r Ry +i 2 ( P r ‘R^ (4,26) 


/BT a 2 -k 2 

a 2 * <->/^ cl+i ^r ] 


(4.27) 


Note that in the last two relations the first terms are very large, of 
order R^ . Consequently, the trigonometric tangents in 0^ and ot^ 
may be written as 


“l *2 

tan y * tan y 


(4.28) 


where exponentially small terms have been neglected. When equations 
(4.25) - (4.28) are substituted into (4.15a,b), one obtains the follow- 
ing equations in a : 


a f tan } % . + 
O l cotan J 2 


1-i ,.* c ~ 2. 

=r (k -a ) 

2kR U 

V 


(4.29a ,b) 


where 


uA ^ 2 


A 


(4.30) 


Higher order terms on the right-hand side of (4.29a,b) are of magnitude 
3/2 

(1/R^) , and have therefore been omitted. As we let the parameter R^ 

go to infinity, the above equations yield the solutions 0^ which perts^n 
to the acoustic modes propagating in an inviscid medium. Hence, the 
zeroth order terms are given by: 


“on * NTT 


(4.31) 


where the integer N is of the form 


N - 2n 


n * 0 , 1 , 2 , . . . 
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for the symmetric mode of order n , described by (4,29a), and 


N « 2n + l 


n » 0 , 1 , 2 , •«« 


for the antisymmetric mode of order n , described by (4.29b). It is then 
a straightforward procedure to derive a more accurate representation of 
Uq by assuming that it may be expanded as: 



Ntr[l+ -—-+£*] 
* 


when N is different from zero, and, 


(4.32) 



(4.33) 


when N is zero. The unknown constants a and b associated with even 
(respectively, odd) values of N , are determined by substituting (4.32) 
and (4.33) into equation (4.29a), (respectively, (4.29b)). The resulting 
expressions are expanded to order 1/R^ « and t * ie coefficient of each 
higher-order term is set equal to zero. The final solutions are then 
given by: 
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(4.35) 


Substitution of the previous results into equations (4.25) - (4.27) yields 
the corresponding expressions for > and /3 ' : 

2 

- + 

(4.36) 
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(k* -(NTT) 2 )(k* +3(Ntt) 2 )^-- ik 4 (-| — +-^) (4.39) 
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The index N may be set equal to zero in (4.36) and (4.37) but not in 
(4.39). When N is an even integer: 

N*2n n * 0 , 1 , 2 , . . . 

the above results describe the symmetric pressure -dominated mode of order 
n , SP(n) , and the associated mode shape is given in Table I. When N 
is an odd integer: 

N * 2n + 1 n-0,1,2,... 

they describe the antisymmetric pressure-dominated mode of order n , 

AP(n) , and the associated mode shape is given in Table II. 

As already explained in Section 3.4, we select in (4.39) the complex 
root which pertains to waves attenuating in the positive z-direction, i.e. 
of negative imaginary part. In the course of the previous derivation we 
tacitly assumed the right-hand side of (4.29a, b) to be much smaller than 
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unity. Furthermore, in order for the above expansions to be valid each 
term has to be smaller than the preceding one. Enforcement of these con- 
ditions results in the following restrictions on the possible values of 
the frequency parameter k : 

« k « R f /3 (4.40) 

R U U 

This double inequality may be written in two other equivalent ways as: 


or 


d 2 \ v * 


4 > l/3 
2 


(4.41) 
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' (JU * 


/ 3 \ 1/2 

(w J y Q ) 


(4,42) 


which is immediately interpreted as an a posteriori justification of the 


name given to the present approximation. Whereas QC ^ is of order unity 
(except for SP(0) ), 0^ and OC are complex numbers of very large 

modulus. Consequently, if we examine the mode shapes of Table I, terms 


such as 


cos « 1 x 
cos (Kj/2 


or 


coaa^x 

cobQL^/2 


may be written in the following manner: 


cos a^x 
cos Qc^/2 


[ot | (x -^) --~|a|(x+j) 

V2 - +e V2 


(4.43) 


1/2 


where |0£j| is very large, of the order of (k R ) 


. Hence, they will 
be of significant magnitude in only a thin layer close to the duct walls 


at x*-l/2 and x»l/2 . In the center of the tube, they will be expo- 


nentially small. As will be clear when we discuss the mode shapes in 


Identical conclusions may be reached with the corresponding expressions 
in Table II. 
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Section 4.6, such • behavior is linked with the diffusion of entropy and 
vortlclty in thermal and viscous boundary layers attached to the walls of 
the tube. If the boundary layer thickness is defined as the distance 
from the duct walls where the magnitude of these terms is e"* of its 
maximum value, the non-dimensional viscous and thermal boundary layer 
thicknesses are respectively: 




kR. L 2kR J 


(4.44) 
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(4.45) 


The lower bound imposed on k in equation (4.40) implies that viscous 

and thermal effects are important only in layers of approximate non- 

dimensional thicknesses and v kP^R* * respectively. By in- 
i' r v 

voking elementary notions of kinetic theory as used in Section 4.2, this 
may be tv*iislated as follows: 


VN 


thN 


« 1 


(4.46) 


i.e., the dimensional boundary layer thickness is approximately equal 
to the geometric mean of the wavelength and mean free path, and it is 
much smaller than the duct width. The high-frequency-wide-tube regime 
might just as well have been called the acoustic boundary layer approx- 
imation. 

To the author’s knowledge, all investiga tions to date with the ex- 
ception of Scarton and Rouleau's have only been concerned with the deter- 
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mination of the complex wave number of pressure -dominated modes, and 
most of them have considered the acoustic boundary layer approximation. 
It is therefore interesting to compare our results with those available 
in the literature. 

Expression (4.38) relative to the zeroth "plane" pressure mode is 
composed of: 

— an inviscid term 

h 

— a second term of order (1/R^) > proportional to the square root 
of the reduced frequency, representing attenuation and disper- 
sion effects brought about by the acoustic boundary layers, and 
first calculated by Kirchhoff (1868) in the case of a circular 
tube . 

— a third term of order 1/R^ > representing higher-order acoustic 
boundary layer attenuation. 

— a fourth term of the same order of magnitude as the preceding 
one, proportional to the square of the reduced frequency, and 
associated with the dilatational attenuation of longitudinal 
waves in the bulk of the fluid. This last part of the complex 
wave number was also derived by Kirchhoff in a study of the pro- 
pagation of plane waves in an unbounded medium. 

Weston (1953a), in an analytical treatment of the characteristics 
of the plane mode in a circular duct, subdivides the wide-tube range 
into two subranges or transition regions. In the wide -narrow tube sub- 
range, k is assumed to be smaller than unity and "closer" to 1/R^ * 
so that terms of the form k/R^ in the brackets of (4.38) can be ne- 
glected . In the wide-very -wide-tube subrange, k is assumed to be 

1/3 

larger than unity and "closer" to R , so that terms of the form 
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or 1/kR^ can be omitted. The simplified expressions re- 
sulting from the manipulation of (4.38) in the above prescribed manner 
are then very analogous to those obtained by Weston in these subranges. 

In order to give the corresponding interpretation of the results 
pertaining to the higher-order P modes, it is convenient to get rid 
of the square root in (4.39) by expansion in powers of 1/R^. Since 
the dominant term in the square root may be positive or negative, one 
must distinguish several ranges of frequencies which are discussed below: 
When k > Ntt , the complex wave number is given by 


- Aw rr 

k -(Ntt) y v 


4i 


(Ntt)' 


k 2 * (k* + (Ntt) 

k -(Ntt) 


2*2 4 

2 k_ k -(Nn) . 1 ik 

2 ' , 2 . 2 ' kR " 2 

k - (Nrr) V 


( y 0 - 1)/ P r + V% 1 ■ 


2 2 
k -(Ntt) 


kR 


V 


(4.47) 


Shaw (1953) studied wave propagation between a pair of infinite paral- 
lei walls and obtained identical results to order (1/R^)^ . In the 
present expression, the expansion has been carried out to higher order. 
This enables us to analyze the different terms in the same manner as in 

(4.38). In particular, the first terra is the well-known inviscid higher 

v 

mode wave number. The total attenuation *^te is of order (1/R^) 2 
and dilatational dissipation is smaller than acoustic boundary layer 
dissipation as can be seen by comparing the second and last term of 
(4.47). In the vicinity of the reduced frequency Nrr , the denominators 
of the perturbation terms in (4.47) tend to zero, so that the above ex- 
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pension becomes invalid. One then evaluates the wave number directly 
from (4.39). 


When the reduced frequency k is equal to the inviscid cut-off 
frequency Ntt , the wave number is: 



1-i 




2NrrR 


p V 0 

1 (— 

V 1 P r 


1 




V 



(4.48) 


In contrast with duct propagation in an inviscid medium, it is not zero. 

Furthermore, its real part is finite, so that a wave still propagates 

along the duct axis. Even though the attenuation rate is now larger, 

of order (1/R ) 4 instead of (1/R , Nrr is not the cut-off fre- 

V V 

quency of the higher-order modes. 

Ntt 

When r < k < Ntt , equation (4.47) is changed into: 
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(Nrr) 2 
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1 

kR 

V 

(4.49) 


It is important to note that the sudden increase of the attenuation rate 
shown in the above relation is not due to enhanced viscous dissipation, 
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but rather to the storage of the acoustic energy in a smaller region of 
space. The physical mechanism is the same as in inviscid wave motion. 
However, the real part of the wave number remains finite and leads to a 
positive phase velocity. 

ft V** 

When k = Ntt or k = y =. — - the real part of the 

[i+(v -i/V ? jy 

U t ft 

complex wave number is almost zero, of order 1/R^ as seen f rom (4.49) . 
Hence, the true cut-off frequency is given by: 



This downward shift in cut-off is solely attributable to thermal effects, 
and disappears when the Prandtl number goes to infinity. The practical 
significance of such a phenomenon is limited, since the attenuation rate 
is already very high at the inviscid cut-off. 

Below Cut-off , the real part of (4.49) leads to a negative phase 
velocity. Wave fronts propagate in the negative z-direction, whereas 
the amplitude is attenuated in the opposite direction. It is shown in 
Appendix C that, in this instance, the acoustic intensity is indeed posi- 
tive, and acoustic energy is moving against the wavefronts. These so- 
called backward-propagating waves were also encountered by Meitzler 
(1965) in the equivalent elastic plate problem, and by Scarton and Rou- 
leau (1973) in a numerical study of the axisymmetric modes in a viscous fluid . 

This completes our discussion of the pressure -dominated solutions 


ft 

In order for the real part of /3 Z to be identically zero to order 1/R^ 
included, one would have to slightly perturb k* around the value Ntt . 
This has not been done here for simplicity. 
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of the dispersion relations (4.15a,b) in the high -frequency-wide -tube 
approximation. 

4,3.2 Entropy-Dominated Modes . Following the methodology de- 
scribed in Section 4.2, we assume the unknown 0^ to be of order 
(1/R^)^ of higher order. In this situation, it is preferable to at- 
tack the problem from the normalized dispersion relations (4.15c,d) 
instead of (4.15a,b). These equations are to be solved for to 


order 


u/y 


3/2 


included* Relations (4.16) - (4.18) enable us to 


approximate Q , OL , and OL by the following expressions: 
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(4.53) 


Since a^ and OL ^ are complex numbers of very large amplitude, their 
trigonometric tangents can be written as: 


d OL 

_ 0 _ 2 
tan — = -tan — = l 


(4.54) 


When these relations are substituted into the dispersion relations 


(4.15c,d), two equations in CL follow: 
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(4.55a,b) 
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^ 1 + 1 (kP rV 
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the upper (respectively, lower) expression pertaining to the symmetric 
(respectively, antisymmetric) entropy -dominated modes. The solutions 
of the above equations are obtained in exactly the same manner as the 
P-modes of the preceding subsection. The final results are given as 
follows: 


“in " Nn 
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(4.56) 


where the integer N is an odd number of the form 


N = 2n + 1 


n = 0 , 1 , 2 , . . 


for the symmetric mode of order n , SS(n) , and an even number of the 
form 

N = 2n n = 1 , 2 , 3 , . . . 


for the antisymmetric mode of order n , AS(n) . Corresponding expres- 
sions for Uq , 0^ , and /S^ are given by: 
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At first sight, * 0 appears to be an exact solution of the anti- 
symmetric dispersion relation (3.71) throughout the entire reduced fre- 
quency domain. However, it is easily checked that the only possible 
mode shape that can be associated with such a solution is identically 
zero. Hence, Oi^ = 0 is not a relevant antisymmetric eigenvalue, and 
as prescribed above, the index n characterizing the antisymmetric 
S-modes takes positive values only. 

The same reasoning as in Subsection 4.3.1 leads to the following 
restrictions on values of the reduced frequency k : 


-L « k « R (4.60) 

R V 

V 


where the upper bound is a consequence of the preliminary assumptions 

made in Section 4.2. Consequently, upon comparison of (4.60) with (4.40), 

it is clear that the range of validity of the S-mode expansions is wider 

than the range of validity of the corresponding P-mode expansions. In 

this investigation, we assume that, by definition, the high-frequency- 

wide-tube approxima tion pertains to the reduced frequency domain 
1 1/3 

— « k « R . The above results are therefore valid, not only in 
R v 

V 


the high-frequency-wide-tube range, but also in the very -high-frequency- 

very -wide -tube range which will be studied in detail in Section 4.5. 

S S 

The coefficients 0^ and a ^ are strikingly similar in form to 
P P 

the coefficients 0^ N and defined in the previous subsection. 

Correspondingly, we will show in Section 4.6 that this mathematical fea- 
ture is linked with a diffusive behavior of pressure and vorticity 
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fluctuations in acoustic boundary layers close to the walls of the tube. 
Expressions for the thickness of these layers may also be derived as in 
Subsection 4.3.1, and they may be shown to be of the same small order 
of magnitude (1/kR^) , as the viscous and thermal boundary layers as- 

sociated with the P-modes. 

The complex wave numbers all have large real and imaginary parts. 

In contrast with the P-modes which have distinct phase velocities and 
attenuation rates, the S -modes have almost identical propagation char- 
acteristics, and are highly attenuated at all frequencies. Nevertheless, 
to each value of the index N , corresponds a specific S-mode shape as 
shown from (4.56). 

Although researchers have been aware of the existence of entropy- 
dominated modes since the earlier work of Chu and Kovasznay, we do not 
know of any previous investigation of their properties for a given 
boundary value problem. In the next subsection, we proceed to determine 
the vorticity -dominated solutions of the dispersion relations (4.15a,b). 

4.3.3 Vorticity -Dominated Modes . The dispersion relations 
(4.15a,b) are to be solved for the unknown 0^ to order 1/R included. 
From (4.16) - (4.18) we may write: 


/3 » -ikR [ 1 - 




V 


kR 


(4.61) 


= 


/ kR . k 2 -**. 2 

a + u/-^[ i-f-s 2 -: 


a, = 


/k(i-p )R ia_ 

( 1 + 1 )/ 2 


2k(l-P r )R y 


(4.62) 


-t 


V 1 


P , R , 

_JL_ ( JL] 
1-P K ? R ' R 

r r 7? V 


(4.63) 
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so that the trigonometric tangents in and 0^ are approximated 

by the following expression: 


tt 0 tt l . 

tan - tan — = i 


(4.64) 


Substitution of the above relations into the dispersion relations 
(4.15a,b) leads to the two following equations in a ^ : 
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OL 

tan •» _2 
cotan J 2 





(4.65a ,b) 


where the upper (respectively, lower) expression corresponds to the sym- 
metric (respectively, antisymmetric) vorticity -dominated modes. By 
making use of the same procedure as in the previous two subsections, we 
arrive at the final results: 

(4 - 66) 

V V 


where the integer N is of the form, 


N = 2n 


n = 1 , 2 , 3 , . . . 


for the symmetric mode of order n , SV(n) , and of the form 


N = 2n + 1 


n = 0 , 1 , 2 , . . . 


for the antisymmetric mode of order n , AV(n) . 

Expressions for a n , a. , and $ are then obtained from (4.61) 

U 1 Z 

(4.63) as follows: 
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The above expansions have the same domain of validity as those describ- 
ing the S-modes, i.e., the reduced frequency k is such that: 


« k « R 


V 


Consequently, these results characterize the vorticity -dominated modes 
in both the high-frequency -wide-tube-, and very-high-frequency-very-wide- 
tube ranges . 

V 

The large real and imaginary parts of the coefficients CKq N and 
V 

a 1N correspond to "diffusion" of pressure and entropy fluctuations in 
thin layers close to the duct walls, the thickness of these layers be- 
ing of the order of (1/kR^) . As seen from a comparison of (4.69) and 

(4.59), the propagation characteristics of entropy-and vorticity-domi- 
nated modes are strictly analogous, and the reader is referred to the 
preceding subsection for further details. 

Scarton and Rouleau (1973) were the first to conduct a numerical 
study of the V-modes (in their terminology, the "B band" of eigenvalues) 
in the case of a circular tube and a non-heat-conducting fluid, and, 
as we shall see in our final discussion, found very similar results. 

They noted that OC^=0 is an exact solution of the axisymmetric disper- 
sion relation, just as it is an exact solution of the symmetric inplane 


dispersion relation (3.70) and numerically derived a corresponding non- 
zero mode shape. However, if a 2 is assumed to be zero in the analysis 
beginning in Section 3.3, one immediately reaches the conclusion that 
the only possible mode shape is zero. Hence, 0^ * 0 cannot be consi- 
dered as a relevant symmetric eigenvalue, and the index N cannot be 
★ 

zero in (4.66). 

We have now completed the determination of the symmetric and anti- 
symmetric solutions of the inplane dispersion relations in the high- 
frequency-wide -tube approximation. At this point, two additional as- 
sumptions which have been implicitly made in the previous derivations 
need to be carefully stated. In order for expressions of the form 
-J kPjRy to be large quantities, the Prandtl number is taken to be 
larger than 1/kR^ . Such a restriction is unimportant since most fluids 
of interest have a Prandtl number of the order of one or larger. Further- 
more, if the S- and V-expressions are to be valid, terms of the form 
Vk(l-P i t )Ry must also be large. Consequently, in these last two 

families, the Prandtl number is in addition assumed to differ from unity 
by a quantity larger than 1/kR^ . This latter restriction could easily 
be removed by examining the particular case where is unity. 

We may now clearly appreciate the advantage of taking , 0^ , 
and a 2 as respective unknowns for the P-, S~, and V -eigenvalues, in- 
stead of the complex wave number. If had been chosen as unknown 

of order or higher order in the dispersion relations, we 

would have obtained the P-eigenvalues only, since they are the only ones 

* k 

A similar situation prevailed in the study of the antisymmetric anti- 
plane modes and antisymmetric S-modes studied in Sections 3.4 and 4.3.2, 
respectively. 
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to be characterized by a wave number of order one. The other two fami- 
lies would have been completely ignored, their wave number being of or- 

Jk 

der as can be seen from (4.59) and (4.69). 

We still have not justified the terminology used to designate the 
different families of eigenfunctions. A definite explanation will be 
given when we examine the characteristics of the mode shapes in Section 
4.6. 


4.4 Low -Frequency -Narrow -Tube Range . 

The solutions in the hlgh-frequency-wlde-tube approximation were 
shown to be valid as long as the frequency parameter k is much larger 
than 1/R^ • For values of k or the order of 1/R^ » the expansions 
break down. In this section, we seek solutions of the dispersion rela- 
tion (4.15a,d) for values of the reduced frequency much smaller than 
1/R^ . The unknown and the dispersion relations are formally expanded 
in powers of k and the range of validity of the results is determined 
a posteriori. 


4,4.1. Pressure and Vorticity-Domlnated Modes . In the high-fre- 
quency-wide -tube approximation, we generally assumed one of the coeffi- 
cients c*q , 0^ , or <* 2 , to be of order unity which lmpliea that the 
other coefficients were very large quantities. In contrast with such 
a situation, we expect in the low-frequency-narrow-tube region all the 
coefficients 0 £q , 0^ , and a ^ to he of the same order of magnitude, 

since terms of the form kR are now very small. Hence, no real advan- 

V 

tage is gained by following the general procedure described in Section 

4.2. In this subsection, we therefore solve the normalized dispersion 

relations (4.15a,b) for the unknown complex wave number /3 and assume 

z 
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/3 to be of order unity. Such a procedure will yield both pressure* 

2 

and vorticity -dominated eigenvalues. The dispersion relations have to 

2 

be expanded to order k inclusive in order to derive accurate solu- 
tions to order k inclusive. From (4.16) - (4.18), a y , 0^ , and a 
may be written in the form: 


a 


0 


10 , [ 1 - 



(4.70) 


a 


a 


1 


2 


ikP R 

is [1+—^] 
1 20 2 


lkR <kR ) 

ifl k [1+ -r2 + — ' ^ 
20 t 2 80 / 


(4.71) 


(4.72) 


When these expansions are substituted into (4.15a,b) and the resulting 

2 

expressions are expanded to order k , terms of zeroth order cancel 
out and the final equations, after division by kR y , and reordering, 
are given by: 
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+ i (^— + ^)(sin i/3 z + i/3 z ) ~ - 0 (4.73a,b) 

r TJ y 

where the upper (respectively, lower) trigonometric function and sign 
correspond to the symmetric (respectively, antisymmetric) dispersion 
relation. From (4.73a,b), we immediately deduce zeroth-order estimates 
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of the eigenvalues of the form 


j8 ■ -i* (4.74) 

2 

where z satisfies one of the following transcendental equations 

sin z « + z (4.75a,b) 

and the + and - signs correspond to symmetric and antisymmetric 

eigenvalues. The solutions of the above equations are investigated in 

detail in Appendix D. They both admit an Infinity of non-zero complex 

solutions. Furthermore, if z is a solution of (4.75a or b) , so is 

★ 

•z and the complex conjugate z . Consequently, if numerical solu- 
tions of these equations are available in the fourth quadrant 
(Re z > 0 , Imz < 0) of the complex z -plane, the solutions in the three 
other quadrants may immediately be obtained by elementary symmetry con- 
siderations . 

Let us Introduce the following notations: the complex root of 

(4.75a) lying in the fourth quadrant of the complex plane and such that 

2nn < Re z < (2n+|)TT n - 1 , 2 , 3 , ... (4.76) 

will be designated z^ . Correspondingly, the complex root of (4.75b) 
which lies in the same quadrant and such that 

(2n+l)n <Rez<(2n + |jTT n « 0 , 1 , 2 , . . . (4.77) 

will be designated z r . The first five roots of each equation have 
been computed numerically, and are displayed in Table III anc' IV. For 
large values of n , they have the following asymptotic forms: 

z + ■ (2n + -biT - i cos h" 1 (2n+~)n as n- + ® (4.78) 

n / l 

z* * (2n +^)n - i cos h* 1 (2n +^)tt as n-+® (4.79) 

nr L 
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In agreement with the "radiation condition" stated in Section 3.4, we 
restrict our attention to complex numbers of negative imaginary part, 
i.e., to solutions of (4.75a,b) lying in the half plane Re z > 0 . 

Hence, the only solutions of interest are {z^ + , ... z^"*",...] and 
{z Q ”,... z and their complex conjugates. The two families, 

{z^ ,... z^ ,...} and { z ^ ,... z^ ,...} have negative imaginary parts 
and therefore correspond to complex wave numbers of negative real parts. 
Since higher-order P-modes were shown to have negative phase velocities 
below their cut-off frequency, in the high-frequency-wide -tube range, 
we naturally associate these two families of roots to symmetric (ex- 
cluding n = 0 1 ,.i antisymmetric P-eigenvalues respectively. Cimilarly, 

since vorticity modes were shown to have positive real parts, we asso- 

+* +* * * 

ciate the conjugate families [z ,... z ,...} and {z ~ ,... z ’ ...} 

in U n 

lying in the first quadrant to symmetric and antisymmetric V-eigenvalues 
★ 

respective! v # 

We do not consider the above reasoning as a rigorous proof of the 
correspondence between high-frequency and low-frequency eigenvalues. 

Such a proof can only be given when the expansions derived in this study 
a* : compared with known numerical solutions of the dispersion relations, 
which allow each eigenvalue to be followed throughout the entire fre- 
quency domain, without any interruption. More specifically, the previ- 
ous identifications will be fully justified when we compare our solu- 
tions with the numerical investigation of Scarton and Rouleau. 

When the zeroth-order or zero-frequency estimates have been deter- 

mined as stated above , it is a straightforward procedure to derive a 
* 

The possibility of identifying "..ese roots with the S-eigenvalues may 
at once be discarded, since such roots still exist when there is no 
heat conduction, where*, the S-modes do not. 
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n 

a 4 * 

n 


i 

7.497676 

2.768678 

2 

13.899960 

3.352210 

3 

20.238518 

3.716768 

4 

26.554547 

3.983142 

5 

32.859741 

4.193251 


TABLE III: 

4 - 4* 4* 

Solutions z + - a. -ip of sin z = z 
n n n 


n 

a “ 

n 


0 

4.212392 

2.250729 

1 

10.712537 

3.103149 

2 

17.073365 

3.551087 

3 

23.398355 

3.858809 

4 

29.708120 

4.093705 


TABLE IV: 

Solutions z = a - iS of sin z = -z 
n n n 
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3?vrf«p«TT 


more accurate representation of the complex wave number by assuming 
that it may be expanded as: 


j8 = -iz(l+ak) 
z 


(4.80) 


where a is determined by substituting (4.80) into (4.73a,b), expanding 
the resulting equations to order k , and setting the coefficient of k 
equal to zero. The symmetric and antisymmetric pressure -dominated wave 
numbers are then respectively given by: 


asp • + r i i lkEl/ k i 

**zn ■ ” lz n [l+ 7« + — - 

4z . 2 n V 

n sin — 


(4.81) 


» ^ ikR iy. , 

^ - -iz-[i+— * + 

zn n , -2 z- R 

4z 2 n y 

n cos — 


(4.82) 


where z and z are the non -zero complex solutions of 


sin z = + z 


which lie in the fourth quadrant of the complex z -plane. 

The symmetric and antisymmetric vorticity-dominated modes are also 
given by the above expansions, where and have been replaced 
by their complex conjugates. Substitution of the results for into 
expressions (4.70) - (4.72) yields the following expansions for a , 


a s? = z + [l +— $ + 
On n , +2 


ikR ./ . j*. , 

4z +2 , l 1 ! '» 

n sin — 


(4.83) 
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(4.84) 


“in ' z n [1+ . +2 


ik(l-2P )R iv, 


r V 


0 k 


— 1 


4z 


. z+ R 
. 2 n y 
sin — 


and 


ikR 


iy r 


. SP + r , y A 0 k i 

a 2n" 'n Cl --« + — rrrr 1 

4z n .in 2 f V 


AP - ikR . ^^0 k 

o£ P * z [1+ % + 0 k 

On n . -2 

4z 


i] 


n cos 


0 2 - R 
2 _n y 


(4.85) 


(4.86) 


AP - (1 - 2P r )R ., iV 0 k . 

“in * Z n Cl + ik + TT R~ ^ 

ln n 4z~ 2 2 z n R y 

n cos — 


(4.87) 



z 

n 


Ci- 


ikR 


4z 


J/ 

-2 




cos 


2 Z n 



(4.88) 


The corresponding V-expressions are obtained by replacing and z r 

by their complex conjugates. 

To derive expansions for the higher-order P-modes and all the V- 
modes, we have assumed the complex wave number to be of order unity. 

In order to determine the characteristics of the only remaining unknown 
eigenvalue, i.e., the zeroth-order SP-mode, is taken to be of the 

form 

0 2 = B k (4.89) 

z 

where the unknown constant B of order unity is found by expanding the 

2 

symmetric dispersion relation (4.15a) to order k inclusive. Terms in 

2 

k cancel out and when the term in k is set equal to zero, we obtain 
the following equation for B : 
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12 iy n 


B= - • 


(4.90) 


V 


so that the SP(0) mode is given by: 





^0- (1 - i) V-R 2: (4 * 91) ; “oo = (1+i) Tr (4 - 92) 


V 



V 


< - 


/kP R 


i^ y q sp 

(l- — 4)(4.93); = (1-iX 

PR 
r v 


r kR 12V 


(1 4 ( 4 . 94 ) 


Upon examination of the ratio of successive terms in the previous 
expansions, one easily determines the range of validity of the low- 
frequency -na rrow -tube approximation : 

k « ~ (4.95) 

V 

As mentioned at the beginning of this subsection, the coefficients 
Gq , , and are now of the same order of magnitude, unity or 

smaller. The associated mode shapes will therefore present smooth vari- 
ations in the cross-section of the tube, and will not exhibit a boundary- 
layer-like character as in the high-frequency -wide -tube regime. 

In the limit of zero reduced frequency, i.e., in the case of steady 
small -amplitude viscous and heat-conducting flow, V-modes and higher- 
order P-raodes are characterized by constant wave numbers. To each 
higher-order P-eigenvalue of a given attenuation rate and negative 
phase velocity, one can associate a V-eigenvalue of identical attenua- 
tion rate and opposite phase velocity. Scarton and Rouleau have de- 
rived similar zero -frequency eigenvalues for the circular geometry and 
Fitzgerald (1972) used the corr sponding eigenfunctions to model the 
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plasma motions in the gaps between successive red cells in narrow-capil- 
lary blood flow. In the present study, additional frequency 'dependent 
terms have been calculated so as to extend the results to a low-frequency - 
narrow-tube region. 

The SP(0)-mode is the well-known unsteady Poiseuille flow analyzed 
by Rayleigh. In contrast with the high-frequency-wide -tube regime, its 
phase speed now varies as the square root of frequency. In the limit 
of zero reduced frequency, this eigenvalue becomes identically zero and 
is therefore irrelevant. Further remarks on the physical meaning of 
these results will be given in Section 4.6. 

4.4.2 Entropy -Dominated Modes . In the above formulation, the dis- 
persion relations were taken in the form (4.15a,b) and such an approach 
failed to yield the characteristics of the S-modes. In agreement with 
the procedure adopted in Subsection 4.3.2, we find it more convenient 
to solve the dispersion relations in the form (4.15c,d) for the unknown 
a^ which is assumed to be of order unity. From (4.16) - (4.18), we 
may write: 


2 2 ikP R (y -1) . R , 

j3 z 2 - -a 1 2 [l+ P V~ P r ( ? _ 'R 2)k ] (4 ‘ 96) 

r rj 


a l" a l 


«n = a, C 1 P - ^)k 2 + 


2 (kP R ) 

r V ^ 


0 ' on 2 5rv 2 r P v V' ‘ cry 2 ‘ o „ 4 

**1 **1 77 **1 8tt l (4.97) 


ik(p -l) (y -l)p r {kd-p )R } 2 

a 2 - “it 1 - -tr^r - *V + zir JL ~ 3 


2a, 


8a, 


(4.98) 


As expected, the coefficients 0^ , a^ , and a^ are seen to be 


of identical order, and values of such that, 
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Ntt 


(4.99) 




a 


1 


where 


N ■ 2n + 1 


or 


N - 2n 


n - 0 , 1 , 2 , ... (4.100) 

n « 1 , 2 , 3 ... (4.101) 


satisfy the dispersion relation (4.15c), or (4.15d) respectively, to 

3 

zeroth order in k . Accurate expressions to order k inclusive are 
needed in order to properly calculate the associated eigenfunctions. 
Consequently, the coefficient Of^ is assumed to be of the form: 

0:^ = Ntt + ak^ + bk^ (4.102) 


where a and b are determined by substitution of (4.102) into (4.15c) 

3 

and (4.15d), and expansion of the resulting equations to order k 

2 

Terms of zeroth order and order k cancel out and coefficients of k 
3 

and k must be set equal to zero. The final results are the following: 





(y -1)(P -l)k ikP R 

“in - »" Cl -- 2 — i Cl 

^ /XT \ *• /XT \ ^ 


2 (Ntt) 


(Ntt) 


(2P -l)y -1 R y (P -l)k 

i r u / 1 IA k. , j U r i i 

+ 1 - —1 ( P c? 5T +i ~Fr J] 

r r T) y tv 


(4.103) 


where the integer N is an odd number of the form (4.100) for the sym- 
metric mode of order n , SS(n) , and an even number of the form (4.101) 
for the antisymmetric mode of order n , AS(n) . Corresponding expres- 
sions for a 0 , a 2 , and 0 2 are given by: 


“on " Nn[l + 


ikP R. (y,-l)P R 0 
2 (Ntt) 2 (Ntt) ?? 


.2 (kP R y 
- jL - T + — ^-] (4.104) 
2 (Ntt) 8 (Ntt) 
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L 


a 


2N 


ik(P -1)R <y 0 -l)P R„ 2 Ck<l-p )R } 

Ntr[l+ — 5 —^ (l-r^)k 2 + V -r JL -l (4.105) 

2(Nn; 2 (Ntt) R rj 8(Ntt) 


£ 


zN 


ikP R (y„-l)P R „ (kP R ) 
•intt c i+— _ i (i--^)k 2 


? ? v R /,v ' a 

2(Nti) L 2(Ntt) *>} 8(Ntt) 


] 


(4.106) 


The range; of validity of these expansions is the same as for the pres- 
sure- and vorticity -dominated modes of the preceding subsection, i.e., 
the reduced frequency k is such that 


k« r 

V 


(4.107) 


In the limit of steady motions, the entropy-dominated eigenvalues are 
purely imaginary and correspond to non-propagating modes of attenuation 
rate Ntt . We do not know of any previous study of these modes, with 
which our present results could be compared. 

We now proceed to an investigation of the dispersion relations in 
the very -high-frequency-very -wide-tube approximation. 


4.5 Very -High -Frequency -Very Wide -Tube Range . 

In Section 4.3 we pointed out that the expansions pertaining to 

the P modes were valid for values of the frequency parameter k , 

1/3 

lower than . Such a requirement is particularly evident if we 

examine relations (4.29a,b) where the right-hand side has been assumed 

1/3 

smaller than unity. For values of k of the order of , the 

high-frequency-wide -tube results break down. Physically, dilatational 

2 

attenuation in the bulk of the fluid, and proportional to k , becomes 
of the same order of magnitude as acoustic boundary layer attenuation 

proportional to . We now reexamine the pressure -dominated modes in 

the domain: 


84 


■r- ^ 


i 


i. 

L 

i 


R^ 1 ^ 3 « k « R^ (4.108) 

where the upper bound on k is a consequence of the continuum hypothesis 
stated in Section 4.2. As mentioned earlier, high-frequency -wide-tube 
results associated with the S- and V-modes are still valid in the pre- 
sent approximation. Consequently, we do not need to consider these 
modes in this section. 

Following a well-established procedure, we solve the dispersion 

relations (4.15a,b) for the unknown which is assumed to be of order 

2 

unity. Equations (4.16) - (4.18) yield expansions of j3^ , 0^ , and 

Oi ^ given by: 


a 2 ,2 r . ..V 1 , V k “o , 
K = k [1 - i( — + F ) 5“-^T ] 

r TJ v k 


(4.109) 


& l = (1-i). 


r kP R . -2 

-E-K [ 1+ i ( ._0_ 
2 L 2 K ? 




R , 

JL 

R V ] 


(4.110) 



(4.111) 


The above expressions are all accurate up to order k/R^ * In (4.109), 

2 

we have included a higher order term in 1/k so as to establish a one- 

to-one relationship between each P-mode and its complex wave number 

0 . However, such a term is negligible in the current determination 

z 

of the equations for . After substitution into the dispersion rela- 
tions (4,15a,b) and appropriate expansion, one obtains the following 
relations : 
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1+i 


(4. 112a, b) 


Of / R 

fcotani _0 _ t 1 + 1 , l kL% 

1 tan j 2 v -1 3l “o 

1 +— — \ K / 

V^r 


Note that the corresponding high-frequency-wide-tube equations (4.29a,b) 

reduce to (4. 112a, b) when the term in 1 /-yj kR^ is neglected in com- 

3 3s 

parison with the term in (k /R^) 3 on the right-hand side of (4.29a,b). 
From (4. 112a, b), it is straightforward to solve for 0f Q . The final 


results are given by: 


tt 0N “ NTT[l+ ‘ 


i« ( 5(f 


1 + 


y ±lWJ 


] 


(4.113) 


V7 


where the integer N is an odd number 


N * 2n - 1 


n = 1,2,3, 


(4.114) 


for the symmetric higher order mode SP(n) , and an even number 


N * 2n 


n * 1,2,3,... 


(4.115) 


for the antisymmetric higher order mode AP(n) . Equations (4.109) - 
(4.111) then yield the expansions of Of^ , CK^ , and : 


r kP k y -2 R 

“in - < 4 - U6) 

r T) v 


tt 2N “ ^ 1 " i ^ 





(4.117) 


2R 


r «* ■ k[l 2 ( F t \\ 2k * <w) V‘ k 2 UV 

i ‘ r-vr 
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(4.118) 


In order to derive the expansions pertaining to the SP(0) and AP(0) 

3 h 

modes, the coefficient 0^ is assumed to be large of o der (k /R r , 
and such that: 


* 


“o ' (1+1)a( » > 

V 


(4.119) 


In this case, equations (4. 112a, b) are shown to still be valid. Fur- 
thermore, since is of the form (4.119), its trigonometric tangent 

is: 

QL 

tan -y* = i (4.120) 


and equations (4. 112a, b) immediately yield the following expressions 
for : 


a 


SP 

00 


AP V 1 k"" 

% ■ (ltl)|1+ 7= )( 2T» 

P v 


h 


(4.121) 


Corresponding expansions for and CC are given by (4.116) and 

(4.117), and the complex wave number is 


« - 1 


z0 


, p - kti-Au+V 1 

zO K L 1 2 


A/ 


k _ i r y °' 1 -f % k 1 
R t, 2 P VV 


(4.122) 


As seen from the second term in (4.118) and the second and third term 
in (4.122), all P -modes are characterized by an attenuation rate pro- 
portional to the square of the reduced frequency, and Kirchhoff's 
free space attenuation term is seen to be a significant part of the 
total attenuation rate. The three coefficients 0 £q , , and Oc ^ at- 

tached to the SP(0) and AP(0) modes are all complex numbers of very 
large imaginary part. Hence, one may consider, in addition to viscous 
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and thermal boundary layer thicknesses associated with 0^ and 0^ 

and defined in equations (4.44) - (4.45), a new layer of thickness of 
, 3 * 

order (R /k ) associated with 0^ , and where pressure fluctuations 
are significant. Such a behavior will be clearly exhibited when we 
examine the mode shapes in the next section. Zeroth order terms simi- 
lar to those of (4.113) were derived by Elco and Hughes (1962) for the 
higher-order P-modes in a viscous fluid contained in a cylindrical duct. 
Weston (1953a) in his investigation of the characteristics of the zeroth 
axlsymmetric sound mode, also finds expressions which are analogous to 
the present SP(0)-mode results. 


4.6 Mode Shapes . 

In the last three sections, the eigenvalues, solutions of the in- 
plane dispersion relations, were analyzed in detail in the low, high, 
and very -high frequency ranges. To each eigenvalue corresponds a given 
mode shape, i.e., given variations of the physical variables with dis- 
tance x from the duct middle plane. Since the symmetric inplane 
eigenfunctions are closely related to the axlsymmetric eigenfunctions 
of a cylindrical duct, we will restrict the study to the determination 
of the symmetric inplane mode shapes. Furthermore, symmetric and anti- 
symmetric eigenvalues were shown to exhibit similar characteristics and 
a discussion of the antisymmetric eigenfunctions would not reveal any 
new physical features. 

The symmetric mode shapes of Table I will be normalized by re- 
quiring that 

p ' (x*0) - 1 (4.123) 

for the SP-modes, 
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S'(x-O) = 1 


(4.124) 


for the SS -modes, and 


d n* 

-j^(x-O) = 1 (4.125) 


for the SV -modes. 

Upon enforcement of these conditions in the relations in Table I, 
the unknown constant A' takes the following forms: 


id 

°o a o k 


( 1 - 


, i 

kP rV COS V 2 




(4.126) 


( 1 - 


kP R ; cos Q£ _ / 2 
TV 1 


in the case of SP-modes, 
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P R 
r y 


<V 1)c p B 2 p 2 

p 0 p l p 0 


(4.127) 


1 


a 2 costt /2 

*1 0 


cos Ofj/2 


in the case of SS -modes, and 


_d 

kR 


cos Qt^ / 2 


V k 


7. 


T 


(4.128) 


V V 


in the case of SV-modes, 

As mentioned in Section 4.1, rhe SP-, SS and SV-mode shapes are 
then obtained in the three frequency and duct-width regir.es Vv substi- 
tuting the expressions for CC^ , 0^ , 0^ , and ^ derived in the pre 
ceding sections into the relations of Table I. 

In order to exhibit in a relatively simple mathematical form the 
main features of each family in each reduced frequency domain, all the 
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terms in the resulting equations can be expanded in powers of the small 
parameter characterizing the frequency range under consideration, so 
as to yield first approximations of the physical variables which are 
uniformly valid throughout the entire cross-section s x * A| . Even 
though we restrict our efforts to the derivation of first-order approx- 
imations, the procedure is long and tedious. Since no mathematical dif- 
ficulties are involved in the derivation, only the final results are 
presented as shown in the next few pages. Tables V, VI, and VII refer 
to the SP-, SS-, and SV-mode shapes in the high-frequency -wide-tube 
range; Tables VIII, IX, and X to the same mode shapes in the low-fre- 
quency-narrow-tube range; and, Tables XI and XII to the SP(0) and higher- 
order SP-mode shapes in the very-high-frequency-very -wide-tube range. 

Alternatively, the complete expansions of oc^ , 0£^ , Of and ^ 
can be inserted into the relations of Table I. The resulting expres- 
sions are very complicated and do not need to be expanded. The mode 
shapes are evaluated directly by computing numerically their variations 
with distance x from the duct axis for given values of the non-dimensional 
parameters, y^ , P^ , , and k . The complex mplitudts of 

pressure, entropy, vorticlty, axial and transversal velocity, and tem- 
perature can then be plotted for the first two modes of each family as 
shown in Figures 5-10. In any figure, each row of rectangular plots is 
attached to a given mode indicated at the left end of the row, and each 
column to a specific physical variable indicated at the bottom of the 
column. On any individual plot, the variations of the real and imagi- 
nary parts of the corresponding physical variable are represented ver- 
sus transversal distance x (vertical axis on the plot) from x = 0 to 
x ■ , i.e., in the upper half of the duct cross-section only. The 
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Table V. High-frequency wide-tube range. SP-mode shapes 
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Table VII. High -frequency -wide -tube range. SV-mode shapes 







cos z 
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cos(2n+l)rrx 



narrow tube range. SS-mode shapes 





cos z 
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Table X. Low-frequency narrow-tube range. SV-mode shapes 











Table XI. Very-high-frequency very -wide -tube range. SP(0)-mode shape. 
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Table XII. Very high frequency very wide tube range 
SP(n) mode shapes, n>0 . 
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Figure 6. High-Frequency Wide-Tube Range. First-order symmetric inplane mode shapes. 
Same values of the parameters as in Figure 5. 
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First-order symmetric inplane mode shape 
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Figure 10. Very-High-Frequency Very -Wide -Tube Range. First-order SP-mode shapes. 
Same values of the parameters as in Figure 9. 














lower half of the mode shape would be obtained by symmetry with respect 
to the duct middle plane for all physical variables except vorticity 
and transversal velocity which are antisymmetric. In the following sub- 
sections, we successively discuss the main results concerning the SP-, 
SS-, and SV-mode shapes in the low, high, and very-high frequency re- 
gimes as they are summarized in the aforementioned tables and graphs. 

4.6.1 High -Frequency -Wide -Tube Mode Shapes . As it is clear from 
Table V, the SP-mode shapes in the central region of the duct are, to 
a first approximation, the same as the corresponding inviscid ones. 
Enforcement of the boundary conditions at the walls is responsible for 
the presence of diffusive terms in the form of decaying complex exponen- 
tials. Such terms are negligible outside thin viscous and thermal 
boundary layers of respective thickness (2/kR^)^ and (2/kP^R^)^ f 
but they play a critical role inside these layers as exemplified in 
Figures 5 and 6. The viscous boundary layer is associated with dif- 
fusion of vorticity of order (k/R^) whereas the thermal boundary 
layer is characterized by diffusion of entropy fluctuations of order 
unity. We remarked in Section 3.4 that a standing wave form, cos 2nnx 
travelling along the positive z-direction, could be viewed as the super- 
position of two travelling plane waves propagating in symmetric direc- 
tions with respect to the duct axis and reflecting against the bound- 

SP 

aries. Let us isolate in the expression for p^ (x) , the plane 
pressure wave 

i(t-nnx-j3 z) 

p'~e Z (4.129) 

As seen from Table V, this wave, upon reflection at the upper wall, will 
give rise to two secondary waves: a highly attenuated "vorticity wave" 
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of the form 



e 



(7 - x) 



(4.130) 


propagating in a direction almost perpendicular to the plane of the 
boundary, and associated with axial velocity fluctuations of the same 
form, and a highly attenuated "entropy wave" given by 


S' ~ e 




kP R 
r y 


(•| - x) 


i[ t+ V 


kP R 
r v 


x - 0 z ] 
z 


(4.131) 


also propagating in a direction nearly perpendicular to the vail, and 
associated with temperature and density fluctuations of the form (4.131). 
The same interpretation may be given for the plane wave propagating to- 
wards the lower boundary. Hence, we may conclude that an SP-mode is 
nothing but the superposition of two plane pressure waves which propa- 
gate in symmetric directions with respect to the middle plane, and upon 
reflection at the duct walls, give rise to very attenuated secondary 
vorticity and entropy waves. 

We also note that SP(0) has a non-zero transversal velocity 
of order (k/R^)^ which varies linearly with x in the central region 
of the duct and rapidly drops to zero in the boundary layers under the 
cancelling effect of viscous and thermal diffusive terms. In contrast 
with this situation, the transversal velocity of the higher-order SP- 
modes is, to a first approximation, the same as in the inviscid case. 

As mentioned in the introduction to this section, the graphs con- 
stitute a more accurate representation of the eigenfunctions since the 

contributions of higher order terms in a , OL , Cr , and $ are 

U 1 Z z 

taken into account. In order to illustrate this feature, we have 
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displayed in Figure 5, the non-zero imaginary part of the SP(0) pres- 
sure amplitude which, according to Table V, is to a first approximation 
identically zero. From a comparison of Figures S and 6, and Table V, 
we infer that higher-order terms only bring about slight changes to the 
mode shapes. 

If we now turn our attention to the results of Table VI and the 
graphs of Figures 5 and 6 pertaining to the SS -modes, we clearly see 
that the perturbations in the central core of the duct are composed of 
entropy waves of order unity associated with temperature and density 
fluctuations of the same order, pressure fluctuations of order k/R^ > 
and velocity fluctuations V' and V* of order (k/R and 1/R , 

respectively. As in the case of the SP modes, enforcement of the 

boundary conditions leads to the diffusion of vorticity fluctuations 

k 2 % 

of order (k/P R ) in layers of thickness ( w~jTp Yr ~ ) close to 

x V ' " r v 

the duct walls. In addition, these entropy waves give rise to diffu- 

i 

sion of pressure fluctuations in layers of thickness (2/kP r R J 2 . 

Such a phenomenon is not directly apparent in Table VI, except through 
the associated transversal velocity fluctuations, because it is of 

SS 

higher order than the term shown in the expression relative to p^ 
on Table VI. However, it definitely occurs, as seen from the more ac- 
curate variations of pressure amplitude displayed in Figures 5 and 6. 
The results may be interpreted in exactly the same manner as those per- 
taining to the pressure -dominated eigenfunctions: An SS-mode can be 

considered as the superposition of two plane entropy waves which pro- 
pagate in symmetric directions with respect to the duct middle plane, 
and upon reflection at the walls of the tube, give rise to secondary 
pressure and vorticity waves. 
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The vorcicity -dominated mode shapes described in Figures 5 and 6 
and Table VII, can be analyzed in a strictly analogous fashion. In this 

case, a vorticity wave associated with axial and transversal velocity 

k 

fluctuations of order 1/kR^ and (1/kR^) , respectively, propagates 

along the axis, in the central region of the tube. In order to satisfy 

the boundary conditions, such a wave has to be supplemented by two types 

of boundary layers. On the one hand, entrop} fluctuations of order 
3/2 

(1/kRy) associated with density and temperature fluctuations of the 
same order of magnitude, are diffused away from the walls in layers of 

thickness (2/k(l-P^)R^)^ . On the other hand, diffusion of pressure 

i c 

fluctuations of order (1/kR^) occurs in layers of thickness (2/kR^)' 2 . 

These fluctuations are also associated with temperature, density, and 
velocity perturbations as seen in Table VII. The SV-modes can also be 
interpreted in terms of two primary vorticity waves and secondary pres- 
sure and entropy waves. 

As seen rrom this detailed discussion, there is ample evidence that 
the high-frequency-wide-tube range is nothing but the so-called acoustic 
boundary layer approximation. The peak of the SP(0) axial velocity in 
the boundary layers, or so-called ’'Richardson annular effect, 11 is seen 
to be the particular manifestation of a much more general phenomenon 
caused by the secondary decaying waves generated at the walls of the 
tube. Such a peak occurs for all higher-order SP-modes as well as for 
the SV- and SS-modes. The amplitude of the SP- and SV -temperature fluc- 
tuations also presents a maximum near the walls. Furthermore, in light 
of the interpretation of the results given above, we have fully justi- 

•k 

cf Richardson and Tyler (1929). 
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fled the terminology used in designating the three families of modes. 

4.6.2 Low-Frequency-Narrow-Tube Mode Shapes . It is clear from 
the expressions of Tables VIII -X and the graphs of Figures 7 and 8, 
that the low-frequency-narrow-tube mode shapes are drastically different 
from their high-frequency -wide -tube counterparts. The amplitude varia- 
tions of the physical variables are now smooth throughout the cross- 
section and no boundary layers can be isolated near the duct walls. 

The vorticity and entropy fluctuations of the SP-modes, for instance, 
can no longer be restricted to thin diffusive layers. They are spread 
throughout the tube and are of the same order of magnitude as the pres- 
sure fluctuations. Similar remarks can be made about the vorticity- 
dominated eigenfunctions. This feature would seem to invalidate the 
argument developed in the preceding subsection regarding the terminology 
used in this study. However, since most practical cases fall into the 
high-frequency -wide-tube range, where the three families assume clearly 
distinct features, we find it convenient to stay with the present con- 
vention. 

The SP- and SV-mode shapes are seen to be closely related, since 
their shape factors z+ are complex conjugates of each other. Further- 
more, on account of the close proximity of the walls, both families are 
characterized by weak temperature fluctuations of order k , so that 
they may be considered as propagating almost isothermally. The para- 
bolic axial velocity profile of the SP(0)-mode is easily identified as 
pertaining to the classical unsteady Poiseuille flow situation. 

The SS-mode shapes are characterized by entropy, temperature, and 
density fluctuations of order unity. The variations of entropy, in par- 
ticular, are unchanged when compared with the corresponding high-frequency- 




wide-tube results. The transversal velocity of order k has not been 

4 

determined, as it would require the knowledge of terms of order k in 
the expansion of 0 t . 

We may therefore conclude that, in the low-frequency-narrow-tube 
range, viscous and thermal effects are dominant in the entire cross- 
section of the tube. 

4.6.3 Very-High-Frequency-Very -Wide-Tube Mode Shapes . As shown 
in Table XI and XII and in Figures 9 and 10, the SP-mode shapes give 

rise to entropy and vorticity diffusion in viscous and thermal layers 

k 

of the same thickness (2/kR^) ' as in the high-frequency approximation. 
However, the SP(0)-mode is no longer a plane pressure wave, as it was 
in the other frequency regimes. More specifically, pressure fluctua- 
tions as well as all other physical quantities are exponentially large 

3 ^ 

in layers of approximate thickness (2R ) close to the boundaries 
of the tube. Outside these layers, a more precise evaluation would 
indicate that the pressure amplitude is of order unity. We also note 
that such layers are thicker than their viscous and thermal counter- 
parts. At very high reduced frequencies, all fluctuations of signifi- 
cant magnitude are concentrated near the walls, resulting in an "annu- 
lar" SP(0) eigenfunction. 

The higher-order SP-mode shapes are pressure released, i.e., the 
pressure amplitude is to a first approximation, equal to zero at the 
walls, as seen from Table XII. Consequently, in contrast with the high- 
frequency case, fluctuations in temperature and axial velocity naturally 
decrease to zero at the boundaries whereas transversal velocities reach 
a maximum close to the walls, and require the generation of secondary- 
entropy and vorticity waves. The pressure waves are associated with 
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entropy fluctuations of order k/R^ in the core of the tube, which may 
be interpreted as increased dissipation in the bulk of the fluid, due 
to the very high frequencies under consideration. The SS- and SV-mode 
shapes do not have to be discussed since their characteristics are the 
same as in the high-frequency -wide -tube regime. 


4.7 Dispersion and Attenuation Characteristics . 

The expansions of the complex wave number pertaining to the SP-, 

SS-, and SV-modes were discussed in detail in Sections 4.3, 4.4, and 

4.5 for high, low, and very-high values of the frequency parameter k . 

From these expressions, one can determine the dimensionless attenuation 

rate, i.e., the attenuation rate per unit diameter along the duct axis, 

| Im /3 | , as well as the phase velocity non-dimens ionalized with respect 
z 

to the speed of sound: 


V = — 

Ph RejS z 


(4.132) 


Both of these parameters are characteristic of the dispersion and atten- 
uation properties of each mode, and have been plotted versus reduced 
frequency k in Figures 11 and 12 for given values of the parameters 


y. , P , R , and R 
0 x V T) 

have been represented. 


The first three symmetric modes of each family 


1 


All curves break down around k “* — , on account of the limited 

V 


range of validity of the low and high frequency expansions. The SP- 

1/3 

results present an additional breakdown around the value k ~ R y 
which separates the high-frequency region from the very-high-frequency 
region. 

In the low-frequency -narrow-tube range, all modes with the exception 
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Figure 12. Symmetric Inplene Modes. Phase velocity versus 

reduced frequency k . Same values of the parameters 
as In Figure 11. 



of the unsteady Poiseullle mode SP(0) , have high attenuation rates in- 
dependent of reduced frequency. As mentioned in Subsection 4.4.1, 

higher order SP- and SV-modes are characterized by almost equal attenua- 

* 


tion rates and opposi'-e phase velocities. Entropy- dominated eigen- 
functions are non-dispersive as their decay rate and phase velocity is 
independent of k . The SP(0)-mode, on the contrary, is highly dis- 
persive, both its attenuation and phase velocity varying as 

Basic changes in the frequency dependence of these two parameters 
take place as we examine the high-frequency-wide-tube range. The at- 
tenuation rate of all SV- and SS-modes increases with-y/k”, instead of 
being constant, and merges asymptotically into the SS(0) and SV(1) at- 
tenuation rates respectively. The same behavior holds for their phase 
velocities. As expected, SP(0) wave fronts propagate with a velocity 
approaching the isentropic speed of sound a^ , and decay at a rate pro- 
portional to -^"k . The most interesting features that occur in this 
frequency range are associated with the higher-order SP-modes. They 
are highly attenuated waves of negative phase velocity below their 
cut-off frequency and weakly attenuated waves of positive phase velo- 
city above cut-off. As the reduced frequency k increases through the 
1/3 

value , the decay rate changes from a ^/k frequency dependence 

2 

due to dissipation in the acoustic boundary layers, to a k depen- 
dence due to dissipation in the bulk of the fluid. Finally, we note 
that the phase velocity is infinite at cut-off, as a standing wave 
pattern is formed in the transversal x-direction. 

It is interesting to compare the results of the present investiga- 


^ote that the absolute value of the phase velocity is plotted in Figure 12. 
This phenomenon was investigated in detail in Subsection 4.3.1 and in 
Appendix C. 
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Cion with those pertainlag to the propagation of sound In an Invslcid 
fluid contained In a two-dimensional duct. As we let the parameters 
Ry , , and go to Infinity, the SP- and SV-famllles vanish alto- 

gether, as their attenuation rate becomes Infinite. The SP-modes then 
coincide with the symmetric sound modes encountered In classical duct 
acoustics, their mode shape being given by 

qp 

p' (x) » cos2nirx (4.133) 

n 


The zeroth mode propagates unattenuated with the lsentroplc speed of 
sound a Q . Higher-order modes are non-propagating below their cut- 
off frequency 2nTT and their decay rate Is 



-\J 4(nrr) 2 -k 2 


(4.134) 


Above cut-off they propagate with no attenuation and their phase velo- 
city is 



k 



2 2 

k -4(nnr 


(4.135) 


Aa we have seen in the high-frequency -wide -tube regime, the presence of 
viscosity and heat conduction results in a finite negative phase velo- 
city below cut-off and a small but non-zero attenuation rate above cut- 
off. It also leads to a dispersive SP(0)-mode with a finite decay rate. 
Furthermore, the SP-characteristics were shown to be drastically altered 
aa we investigated the low- and very-high-frequency regimes. In these 
two ranges, the invlscld results do not lead to meaningful predictions. 


4.8 Comparison With a Numerical Study . 

It la of Interest to compare the analytical results derived in this 
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chapter with those obtained by the numerical approach of Scarton and 
Rouleau (1973). These investigators assume that the fluid is viscous 
but non-heat-conductlng, and that it is contained in a rigid cylindri- 
cal duct of circular cross-section. Their work is concerned with the 
determination of the characteristics of the axisymmetric modes by the 
method of eigenvalleys. Since it was pointed out in Section 4.2 that 
the dispersion relations pertaining to the axisymmetric and lnplane 
symmetric eigenvalues respectively, could be written in analogous 
mathematical form, it is legitimate to seek the common features be- 
tween the numerical solutions of the first relation and the analyti- 
cal solutions of the second relation. 

For convenience, we discuss Scarton and Rouleau's results in terms 
of the three non-dimensional parameters k , R^ , and R^ introduced 

1 , 

in Section 3.2, where the duct width d is replaced by the radius R . 
As we let the Prandtl number go to infinity, the attenuation rate of 
the entropy -dominated modes becomes infinite in both geometries so that 
the SS-modes and their axisymmetric counterparts are identically zero. 
Hence, this section is concerned with a comparative evaluation of only 
two families, namely, the pressure- and vorticity-dominated modes. 

These two families are referred to by Scarton and Rouleau as the A- and 
B-bands, respectively. In Figures 13 and 14, we have reproduced typical 
plots of their dimensionless attenuation rate, l.e., attenuation rate 
per unit radius, and dimensionless phase velocity versus reduced fre- 
quency k , for a given set of values of the parameters and R^ . 

★ 

In their formulation, Scarton and Rouleau prefer to use the equivalent 

set k , 1/R^ , and Kg/p^R . 
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Absolute value of dimensionless spstisl attenuation factor, \y f \ 



Figure 13. (From Scerton end Rouleeu, Fig. 4(b), p. 604) 
Dlmenslonlesa attenuation rate versus reduced 
frequency k-tuR/a* . Non-heat-conducting vis 
cous fluid; rigid circular duct. R “a-R/v ■ 
10 , K q - 0 . U 
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Absolute vftluo of dim«M>on)ea« phMe velocity, \<*\ 



Dimensionless frequency, /’ 


Figure 14. (From Scarton and Rouleau, Fig. 5(b), p. 608). 

Dimenaionleaa phase velocity versus reduced fre- 
quency k . Same values of the parameters as in 
Figure 13. 
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Figure 15. (From Scarton and Rouleau, Fig. 7, p. 613). 

Zeroth and first order mode shapes for k = 10 , 

R =10^ . 0 , real part; A imaginary part. 
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The reduced frequencies considered in these grap' s span the entire high- 
frequency range and the edges of the low- and very-high-frequency re- 
gions. The overall similarity between Figures 13-14 and Figures 11-12 
is striking. In the low-frequency regime, the higher order A- and B- 
modes have identical attenuation rates and opposite phase velocities, 
in complete agreement with the observations of Section 4.7. As the re- 
duced frequency k is increased through the transition region, the 
propagation characteristics of the A- and B-band closely approach those 
of the SP- and SV-modes in the high frequency range. This observation 
confirms that in the low-frequency region, the solutions of equation 
(4.75a) located in the first and fourth quadrant of the complex plane 
do correspond to the SV- and SP-eigenvalues, respectively. As men- 
tioned in Section 4.4, such a choice could not have been fully justi- 
fied without a knowledge of the behavior of the SP- and SV-complex wave 

numbers in the transition region between low ana high frequencies. In 
the high-frequency region, higher-order A-modes experience in the same 
manner as higher order SP-modes, a sharp decrease in decay rate and a 
change of sign in phase velocity, as k crosses their cut-off frequency. 
Since Scarton and Rouleau consider a purely viscous fluid, the cut-off 
frequencies are the same as in the case of an inviscid medium, and no 
downward shift is observed. As the reduced frequency is further in- 
creased, the yr dependence of the SP-attenuation rates gradually 
2 

shifts to a k dependence. However, there seems to be an apparent 
discrepancy between the respective magnitudes of the AO- and SP(0) - 
decay rates above the cut-off frequency of the first-order pressure- 
dominated mode. Whereas the axisymmetric AO-mode is more attenuated 
than the higher-order A-modes in this frequency domain, the two-dimen- 
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sional SP(0)-mode is less attenuated than the higher-order SP-modes. 
Beatty (1950) used an equivalent impedance model to show that the 
higher-order pressure modes were indeed the least attenuated of all 
axisymmetric modes, and that the opposite situation prevailed in the 
two-dimensional case. A study of the cylindrical configuration on the 
same lines as the present investigation would confirm this particular 
point. 

In Figure 15, we have reproduced plots of the modal amplitudes per- 
taining to the zeroth and first modes in the A- and B-band, as obtained 
numerically by Scarton and Rouleau. The values of k and corre- 

spond to a typical high-frequency-wide -tube situation. Since the real 
and imaginary parts have been normalized by their maximum in the inter- 
val O^R^l , a detailed comparison is impossible. However, we note 
the general similarity between Figure 15 and Figures 5 and 6 of the pre- 
sent study. In both instances, pressure- and vorticity -dominated modes 
are characterized by the diffusion of vorticity and pressure, respec- 
tively, in thin boundary layers close to the walls. Scarton and Rouleau 
note that the thickness of these layers increases as the mode index in- 
creases and the frequency parameter decreases. This trend is clearly 
demonstrated by the analytical expressions (4.44) and (4.45). 

Even though the value of the reduced frequency associated with the 
mode shapes of Figure 16 is located at the edge of the low-frequency- 
narrow-tube region, we note in these plots the absence of boundary layers. 
Pressure and vorticity fluctuations diffuse in the entire cross-section, 
a situation which is very reminiscent of the low-frequency mode shapes 
displayed in Figures 7 and 8. 

As is clear from the examination of Figures 13-16, Scarton and 
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Rouleau claim that there exists a BO-eigenvalue associated to an eigen- 
function which is not identically zero. If this were to be true, we 
would have to admit that the corresponding exact solution of 

the inplane symmetric dispersion relation (3.70) is a relevant eigen- 
value with a non-zero SV(0) -eigenfunction. As mentioned in Section 
4.3.3, if a 2 is assumed to be zero in the basic equation (3.42), one 
is naturally led to conclude that the only possible mode shape is iden- 
tically zero. Hence, the solution 0^ = 0 and its axisymmetric counter- 
part BO are not genuine eigenvalues. 

4.9 Concluding Remarks . 

The discussion of the last three sections was restricted to the 
symmetric modes. The propagation characteristics and shapes of the 
antisymmetric modes could be analyzed in an identical fashion, in light 
of the solutions of the antisymmetric dispersion relation derived in 
Sections 4.3 - 4.5. 

The determination of the inplane modes, in view of the complexity 
of the basic relations describing their motion, required the implementa- 
tion of perturbation procedures. We were therefore unable to obtain 
analytical results in the transition regions where no small parameter 
could be defined. No such difficulties were encountered in the case of 
the antiplane modes and their characteris tics were determined exactly 
in the entire reduced frequency domain. A simple explanation of this 
difference between the two types of motion can now be given: the anti- 

plane modes are solely composed of vorticity waves, whereas the inplane 
modes are the result of the coupling of primary waves and secondary 
waves generated at the boundaries. In particular, the inplane vorticity- 


dominated modes, in sharp contrast with their antiplane counterparts, 
give rise to pressure and entropy fluctuations. 

The following question remains to be answered: given a source dis- 

tribution located inside the duct, in the cross-sectional plane z-0 , 
what is the nature of the disturbances propagating along the duct axis. 
Their character will naturally depend on the modal content of the source 
amplitude for each reduced frequency k . A generator of pressure fluc- 
tuations, such as a piston oscillating in the axial direction, or a 
pulsating sphere, will give rise to a wave pattern composed of pressure - 
dominated modes. Similarly, if a heat source or a distribution of heat 
sources is situated inside the tube, the wave will consist of entropy- 
dominated modes. Finally, a piston oscillating in its own plane in the 
x- or y-direction will give rise to a wave pattern made up of inplane 
vorticity-dominated modes or antiplane vortical modes. The number of 
eigenfunctions to be considered and their respective weight will be de- 
termined by the modal structure of the source. 

We noted in Chapter II that pressure and entropy waves are inher- 
ently coupled as they propagate in an unbounded medium, whereas vortical 
waves are not. From the results of this chapter, we may conclude that 
the presence of solid boundaries has led to additional coupling between 
vortical fluctuations and the two other types of fluctuations. For in- 
stance, the presence of turbulence at some station along the duct will 
generate vorticity-dominated modes which include pressure and entropy 
fluctuations. In the absence of boundaries, the radiation field would 
solely be composed of purely vortical waves. It is important to remark 
that this conclusion is not in contradiction with Lighthill's theory of 
aerodynamic noise generation. The present formulation does not include 


the pressure waves generated by the non-linear fluctuating Reynolds 
stresses. Hence, turbulence does not give rise to pressure disturbances 
in the linear approximation, when the medium is unbounded. However, 
when there are solid walls, we have shown thst even in the linear ap- 
proximation, pressure fluctuations are induced at the walls and propa- 
gate along the axis at a very high attenuation rate. 

In Figure 17, we have summarized on a k - diagram the ranges 
of the parameters k and R^ in which inplane solutions have been ob- 
tained. As explained in Section 4.2, regions of the k - R^ plane 
located on the left of the straight lines R^=l and k = R y have to 
be excluded on account of the continuum hypothesis. The lines k = 1/R^ and 
k * R separate the low- and high-, the high- and very-high-frequency 

ranges, respectively. In fact, our results are not valid for points 
located near these lines, and a more accurate representation would have 

to separate the different domains by transition bands around the lines 

1/3 

k-l/R and k«R /J . 

V v 
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Figure 17. Domain of validity of the low-, high-, 

and very-high-frequency approximations in 

the R -k plane. 
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V. CONCLUSIONS AND RECOMMENDATIONS 


5.1 Conclusions . 

The important achievements of this study are outlined below. 

1. A new mathematical formulation of the small fluctuating motions 
of a viscous, heat-conducting and compressible fluid has been developed, 
which reduces such problems to the determination of three unknown func- 
tions, namely, the acoustic, thermal, and viscous potentials, satisfy- 
ing three linear partial differential equations. As a consequence , the 

total number of unknown functions has decreased from five "physical" 

— ♦ 

variables, i.e., p r , p' ,S' , T 1 , and V 1 , to three "auxiliary" variables, 

$ , 3> , and A . 

a th 

2. This model has been applied to a detailed investigation of the 
propagation of sound in a two-dimensional infinite duct with rigid and 
isothermal walls. It has been shown that four distinct types of wave 
motions, i.e-, inplane and antiplane waves of symmetric or antisymmetric 
amplitude with respect to the duct middle plane, can be considered in- 
dependently. 

3. The characteristics of the symmetric and antisymmetric anti- 
plane modes have been obtained exactly for arbitrary values of the para- 
meters. They consist of purely vortical velocity fluctuations. 

4. In the case of symmetric or antisymmetric inplane motions, 
three families, namely, pressure-, entropy-, and vorticity -dominated 
modes, respectively, have been distinguished and their characteristics 
have been determined in the low-frequency -narrow- tube , high-frequency- 
wide-tube, and very-high-f requency-very -wide -tube approximations . 

5. The analytical results pertaining to the symmetric inplane modes 
have been found to be in good agreement with those obtained in a numerical 
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study of the axisymmetric modes of a cylindrical duct in the limit of 
zero heat-conduction. 

5.2 Recommendations . 

The general features of the fluctuating field generated by an arbi- 
trary source distribution were discussed qualitatively in the last sec- 
tion of the preceding chapter. A mathematical treatment of this prob- 
lem will require a detailed investigation of the completeness of the 
modes and their orthogonality properties. A proof of their complete- 
ness will ensure that any arbitrary disturbance may always be regarded 
as a linear combination of the eigenfunctions. When the scalar product 
of two modes is defined as the integral of their product over the cross- 
section, it is immediately verified that the modes are not orthogonal. 
However, a systematic inquiry might reveal the existence of more com- 
plicated orthogonality relationships. Such properties would greatly 
simplify the calculation of the coefficients multiplying each eigen- 
function in the linear combination mentioned above. 

It would also be of interest to compare the analytical results of 
the present investigation with those obtained by a numerical study of 
the inplane dispersion relations. The analytical results themselves 
provide excellent first guesses of the roots in a wide range of values 
of the parameters, and the implementation of a standard Newton-Raphson 
procedure would enable us to determine numerically the eigenvalues in 
these ranges as well as in the transition regions which have been ex- 
cluded in the course of the perturbation analysis. 

The methodology developed in the present work and the aforementioned 
suggested extensions then could be applied in a straightforward manner, 
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to the determination of the characteristics of the axlsymmetrlc and 
spinning modes in a circular duct. The main difference would reside 
in the introduction of Bessel functions instead of trigonometric func- 
tions. 

In many practical cases, such as the propagation of sound in the 
atmosphere or through aircraft engine-ducts, velocity and temperature 
gradients are present in the medium. It would be worthwhile to examine 
possible extensions of the mathematical formulation of Chapter II to 
include complex mean flow situations. Such a theory would encompass 
a broad range of problems and provide a general framework to analyze the 
propagation of sound in a viscous and heat-conducting fluid which is 
Inhomogeneous and non -uniformly moving. * 


l 


128 




APPENDIX A 

CALCULATION OF (||) , (|^) , (||) 


In this Appendix we will omit the subscript 0 with the under- 
standing that we are considering the state of the medium at rest. The 
characteristics of the fluid are assumed to be specified in terms of: 

- the temperature T 

- the density p 

- the isentropic speed of sound a A * _/ (t^) 

0 -V dp s 

- the specific heat at constant pressure c^ 

- the ratio of specific heats V 

c^d bT 

The partial derivatives (-rfO , (r-) , (^-) are to be expressed 

p s ° S p 

as a function of the five properties of the fluid mentioned above. In 
the following derivation we make use of the First and Second Laws of 
Thermodynamics without introducing any additional assumption such as a 
particular equation of state. 

From the definition of partial derivatives, we may write: 


dp - a 2 dp + (||) dS 


(A-l) 


ax . <0) dp . (§) ds 

S p 


(A-2) 


Equation (A-2) is then used to express dS as follows: 


75t~ [dT ‘ <S> de ] 


(A -3) 


The First and Second Laws of Thermodynamics reduce to the differential 


form: 
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de « TdS + -^dp 
P 

U-4) 

and when (A-3) la substituted into (A-4) we obtain the following 

relation: 

T ^ 

de ■ W 

< as ) “ ( as> 

P P 

(A-5) 

The coefficient of dl in (A-5) la nothing elae but 

the specific heat 

at constant volume so that 



T 

CV ‘ & 
( *p 

(A-6) 

or 

Ms . X 

W c 

P V 

(A-7) 

which provides an expression for as a function of T and c^ . 

In view of 

r 

this result, equations (A-3) and (A-5) may be rewritten as: 


dS - ^ [dT - (|^) do ] 

s 

(A-8) 


de * c dT + [ -% - c (r~) 3dp 
v 2 v ip . 

P 5 

(A-9) 

Since both 

de and dS are exact differentials, we 

have from (A-8) 

and (A-9) : 


(A-10) 


‘it't • 7 #„ • c ^<s< 0 ) > i 

S P 

(A-ll) 
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Substitution of (A-ll) into (A-10) then results into the following 
simple relation: 


c ,*L) . X (*£) 
C v^p ) . 2 V 

S p p 


(A-12) 


where (|^) is to be expressed in terms of (|^) . When (A-8) is used 

P P 

in (A-l), dp is related to dT and dp through the equation 


dp “ rjr (|g) dT+ const, dp 
P 


(A-13) 


so that 


<t?> ■ f<H> 


(A-14) 


P P 

and, with the help of (A-14), (A-12) becomes: 


,3T. 1 

<55> ' 7<3s> 

S p p 


(A-15) 


St 

The above result provides a first equation between the unknowns (r^-) 

3p q 

and (^jr) ♦ The second equation is obtained by expressing the First 

o 

and Second Laws of Thermodynamics in terms of the enthalpy h : 


dh 


TdS+ ^ 
P 


(A-16) 


Equations (A-l) and (A-2), together with the result (A-7), then lead 
to the following relation between dS , and dp and *.T : 


dS 


a 2 dT - (|J) dp 
3p S 


a 2 T 

C 


. ( «£) (^1) 
V 'Sp', 

P £> 


(A-1T) 


which, after substitution into (A-16) results in: 


•2 . A, (Si, 

c v <8S> p M 


dT + const, dp 


(A-18) 


The coefficient of dT in (A-18) is the specific heat at constant pres- 
sure so that 


c = a 2 T 

S <!><!;> 
v p S 


(A-19) 


which provides a second equation between the unknowns (^— ) and (^) 


C (|f) (It) - (y-i)a 2 T 


V*ya P ' s 


(A-20) 


Equations (A-15) and (A-20) are immediately solved for (— ) and 

dp n 

(^) , and the final results are shown below: 


<!> - 


(y-i)T 


(A-21) 


V s » 


(A-22) 


<g> - ? 

P P 


(A-23) 
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APPENDIX B 


DERIVATION OF SPLITTING THEOREM 


Let us first introduce the following notations: 


L, ($ +*5 , ) = + $ .) - — V 2 (<l> + d> ,)+ £- 

1 a th' dt v a th' p Q v a th' p Q 


(B-l) 


r* 3A n 2 t* 

L 2 (A) " SVV A 


(B-2) 


L , ) = $ - — — 

3 th th Vo 



V. 

c T "0 
P 0 ° 


k T' 


(B-3) 


The problem considered in this Appendix may then be stated in the fol- 
lowing manner: 


* _ k _ /V —+~k - 

Assumptions : Let the set 1 , *7 , A j satisfy the following 


system: 


k _ k 


-+k 


grad [ L^ (f^ +* th ) 1+ curl [L 2 (A )] = 0 


(B-4) 


- 0 


(B-5) 


and V' , p 1 , and S' be given by the relations: 


V' = grad + 'f* 4- cur 1 A 


(B-6) 


V 2„2x* 

-r^- = -p a 7 
fit M 0 0 a 


(B-7) 


Ss’ 

at = a o 



F 0 n 2 A * 

<V 1)T o th 


(B-8) 


with 


-tit 

div A =0 


(B-9) 
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Question : Is it possible to find a new set { , A , A } which 

will yield the same values of the physical variables through the rela- 
tions 


V* = grad (A + A..)+curlA 
a tn 



-Vo* ***. 



with 


div A = 0 


(B-10) 


(B-ll) 


(B-12) 


(B-13) 


and which will satisfy the simpler system of governing equations: 


vw = 0 


L 2 (A) = 0 


W - ° 


(B-14) 

(B-15) 

(B-16) 


Proof: Let us examine the following new set { A , A . , A } : 

a tn 


* . * * p * * 1 

A = A +A - l.(A +A^. )dr - 

a a th J 1 r ** u 


V 1 


A. = 



a th ' vo'V c p 0 v° 


k„T* 


k„T'(B-17) 


(B-18) 


A = A* -J L 2 (A' V )dr (B-19) 

0 

and verify that it complies with the requirements (B-10'> - (B-16). 
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(a) From (B-19) we have: 


divA = divA -J divL^CA )dT 
0 


(B-20) 


and 


l ^ ^ | ^ r\ 

divL 2 (A ) - ^-divA - div V A = -^-divA - t/ Q V divA 


Since (B-9) is satisfied, (B-20) and (B-21) require that 


(B-21) 


divA = 0 


(B —22) 


which proves (B-13). 

(b) Let us now consider 


-» * •k 

grad (<P^ +3^)+ curl A = grad(^ +<P t ^)+curlA 


a th' 
t 


-J { grad curl L^(A )}dr 


and therefore, from (B-6) and (B-4) 


(B-23) 


grad(3> +3) )+curlA = V* 
a tn 


(B-24) 


which proves (B-10). 

(c) Let us calculate 


7 2 $ = V 2 * * + V 2 [<f> * - — 

a a th p A a 


V 1 


o"o V c P 0 T o 0 


k A T'] 


-1 72 vV t+ * t h )dT 


(B-25) 


When we take the divergence of (B-4), we have 


V 2 L (<b *+* *) = 0 
la tn 


(B-26) 


135 


so that (B-26) and (B-5) require that 

v 2 * = v 2 $ * 

a a 

which, together with (B-7) proves (B-ll) . 

(d) Similarly, we may write 



so that, from (B-5) we have 

,2 *th - 72 < 


(B “27 ) 


(B-28) 


(B -29) 


which together with (B-8) proves (B-12). 

(e) From (B-17) - (B-18) we may write: 

t 

VW ' L l (4 a* + *th ) - ii I l l < ® a * + V dT 

0 

rj t 

+ — | V 2 L. ($*+$>*) dr (B-30) 

1 a th 

0 0 

so that with the use of (B-26) 

L l (4> a + *th ) = 0 (B-31) 

which proves (E-14). 

(f) From (B-19) we have 

t t 

L 2 (A) = L 2 (A f ) J L 2 (A*)dT + v 0 J 7 2 L 2 (A“)dr (B- 32 ) 

0 o 

When we take the curl of (B-4) and take into account (B-9), we may write 
V 2 L 2 (A*) = 0 ( B - 33 ) 
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which shows that (B-32) implies: 

L 2 ( A) = 0 (B-34) 

(g) Finally, using (B-19), we have 

L 3 ( *th ) = ° (B-35) 

which proves (B-16). 

The new set { A , , A } therefore meets the requirements 

(B-10) - (B-16). In particular, it is governed by the system of partial 
differential equations (B-14) - (B-16). 
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APPENDIX C ACOUSTIC POWER AND ENERGY TRANSPORT 


VELOCITY OF HIGHER-ORDER SP-MODES 


The radiation condition uaed in the present investigation requires 
that waves be attenuated in the positive z-direction, i.e., that the 
imaginary part of the complex wave number be negative. We showed in 
Section 4.3.1 that such >a condition results in the existence of pres- 
sure-dominated backward -propagating waves of negative phase velocity 
below cut-off. In this Appendix, we prove that in spite of the change 
in sign of the phase velocity through cut-off, the acoustic energy of 
the higher order SP-modes always propagates in the positive z-direction. 

Karamcheti (1974) showed that the classical definition of acoustic 
intensity is not affected by the introduction of viscosity and heat-con- 
duction effects. The acoustic power is therefore taken to be given by 

% 

P = R e f f J p ' v ;*dx} (C-l) 

-h 


where the duct width d has been introduced because x is a non-dimen- 
sionalized coordinate* We are interested in deriving a first approxi - 
mation of the acoustic power radiated by the higher-order SP-modes in 
the high-frequency range. In this context, the acoustic power needs 
to be evaluated to order (1/R^)^ inclusive. The pressure and axial 
velocity amplitudes are given in Table I, After expansion of the coef- 
ficients multiplying the cosine functions, p f and can be written 

as 


p* (x,z,t) 


lp o*o k 


cos CKqX 


. R . cosa.x z) 

o r tj v i 

(C-2) 


138 


(x,z,t) 


< A i* * cosa n x u cos< 

. r 2 ify 

d *z L cosC^*/2 ''O C08< 


Vo- 


where the SP superscript and n subscript have been omitted for con- 


venience. In the preceding expressions $ a n , OL , 0£ o , and 0 are 

U 1 2 Z 

given by equations (4.35), (4.36), (4.37), and (4.39), respectively, 
and A 1 by equation (4.126). In order to evaluate the power, we need 
the values of the following elementary integrals 


J 


h * 

cosO^xcosG^x 


^ cos a^/2co3QC*/2 


dx 


2 5^(a 0 tana Q /2 -a 2 tana 2 /2) 

(C-4) 


“0 -°2 


and 


I 


L ^ 

■ cosa 0 xcosa 0 x 2Im{a 0 tana 0 /2} 

dx = — 


•h 


cosa Q /2 cosO*/2 


(C-5) 


ImOf. 


After substitution of the expansions of and 0£ 2 , these integrals 

can be reduced to 


and 



* 

cos G£qX cos OL 2 x 
cosa Q /2 cos a*/2 


dx 



★ 

COSQ£qXCOS X 

cosa Q /2 cosa*/2 



*2 

k -4(nn) 
4(nn) 2 


2 



(C-6) 


(C-7) 


where terms of order 1/R^ have been neglected. Similar expressions 
can be derived for the other integrals which appear in the acoustic 
power. Integrals of the same form as (C-6) and corresponding to the 
paits ( Ot Q ,l X*) , (a 1 ,a Q *) , (a^.O^*) , and (0^,0^*) are found to be 
of order (1/R^ . When these intermediate results are substituted 
into (C -1) only two terms need to be retained, so t * c the power P 
is given by 
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(C-8) 


P 


^ *•*■*- 



COS Q£qX COS a Q X 

cosa Q /2 cosa*/2 dx 


-I 




•A 1 

cosa Q x cosa 2 x 
cosa Q /2 cos a*/ 2 



2ImjS z 
e z 


The first integral represents the contribution to the acoustic power 
of the irrotational part of the fluctuating motion whereas the second 
integral involves coupling between the irrotational and rotational 
fluctuations in the acoustic boundary layers. When the constant A' 
given by (4.126) is expanded and use is made of (C-6) and (C-7), the 
power finally equals 



(C-9) 


Instead of reasoning in terms of acoustic power, it is convenient to 
normalize P by the acoustic energy per unit length along the duct 
axis. In other words, we prefer to consider the energy transport velo- 
city non-dimensionalized with respect to a Q and defined as follows: 


V 

en 


(C-10) 


where E is the acoustic energy per unit length. Note that in the 
present case of a dissipative medium, the energy transport velocity 
does not necessarily coincide with the group velocity. The first ap- 
proximation of E g in the high-frequency range isidentical to the in- 
viscid acoustic energy. Furthermore, in the case of an inviscid medium, 
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Che acoustic energy is of order unity throughout the entire frequency 
domain. Consequently, in order to obtain a first approximation of , 

only the inviscid value of E g is needed. It is given by 


4 Vo 2 ’* 


h d ^ 

d r p o d (*-» - * 

4 


J p'pf*dx+^lj V'.V* dx 


(C-ll) 




where p' and V 1 are the pressure and velocity fluctuations associ- 
ated with the symmetric acoustic modes propagating in the inviscid 
medium. When we use the results mentioned at the end of Section 4.7, 
we find that E g is given by 


and 


*s / . 2 

“Vo 


, ,2 2Im(9 s 

E 1 

s 2 

p o a o 


when k > 2nn 


(C -12 ) 


when k £ 2nrr 


(C-13) 


From equations (C-9), (C-10), (C-12), and (C-13), and the expansions 

for the complex wave number j3 derived in Subsection 4.3.1, one can 

z 

determine the acoustic power and the energy transport velocity, and 
follow their variations as the frequency parameter k decreases. 

When k > 2nn , the acoustic power is given by 


-L-J 

4 p o*o k V 


2 2 
k -4(nn) 4 


(C - 14) 


and the energy transport velocity is 


v . V k 2 -4(nTT)J 


en 


(C-15) 
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These expressions do not exhibit any new features. They are identical 
to the inviscid results. 


When k - 2nrr , the wave number is given by equation (4.48) and the 


power and energy transport velocity respectively equal 

p i /Vjf . 1 * n - 8 nn (v^) ^ * l ”»' 

' 2 Vo (vfJ 2nnR „ C ° a 8 * 


(C-16) 


en 


( y 0~ l \ ( 1 } tt 

(yr) 2n " R p co “ 8 


(C -17) 


The energy propagates with a velocity and decay rate proportional to 

(1/R • In the inviscid case, both these quantities would be zero 

and remain equal to zero for lower values of k . 

2nrr 


When 2nrr < k < 


tion (4.49) and we have: 

P » f 
4p ( 


(i-Kvq-D/VT) 


^ , the wave number is given by equa- 


A A A(4(nTT) 2 -k 2 ) J. k* 2 -4(n") 2 + 2 j g -2 ^(nn) 2 -k' 

0 a 0 k V kR y 4(nTT) 2 -k 2 


(C-18) 


V 

en 


k 

4(nrr) 2 


/ 2(4 (ntr) 2 -k 2 ) r 

V kR L 

V v 


2 

k -4(nTT) 

2 2 
4(nn) -k 


+ 2 ] 


(C - 1 9 ) 


The first and second terms in the brackets of (C-18) represent the con- 
tributions of the irrotational and rotational components of the fluc- 
tuations, respectively. They are positive and of the same order of 
magnitude. The total power is positive. Equation (C-19) can be inter- 
preted in the same manner. 

When k -=z — p , the appropriate values of P and 

(1+O'o-1)/Vp r )'* 

V gn are obtained by simply replacing the reduced frequency by its value 
in (C-18) and (C-19). It is easily checked that, to this approximation. 
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the power Is associated with the coupling of rotational and lrrotatlonal 
fluctuations in the boundary layers. The purely lrrotatlonal pert of 
the power la of higher order. Hence, at cut-off, the energy propagates 
at a high attenuation rate In thr viscous boundary layers. Below cut- 
off , the rotational component is at.* 1 1 positive and larger in magnitude 
than the negative lrrotatlonal component, as seen from equation (C-18). 
In this frequency range, the rotational part of the energy propagates 
in the positive s-dlrection in the vlucous layers, whereas the lrrota- 
tional part propagates in the opposite direction and is distributed in 
the entire cross-section. The net power is -till positive. 

The results of this discussion are summarised in Figure 18. The 
velocity of energy transport pertaining to the SP(1) mode is plotted 
versus the reduced frequency. It is clear from the graph that the vor- 
tical fluctuations play a crucial role in keeping the velocity of en- 
ergy transport positive. Identical conclusions would result from an 
analogous study of the antisymmetric modes. 
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ENERGY TRANSPORT VELOCITY 


10 



Figure 18. SP(1) -energy transport velocity versus 

reduced frequency k . Same values of 
the parameters as in Figure 12. 
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APPENDIX D 


SOLUTIONS OF SIN Z • ± Z 

In this Appendix, we discuss the mein properties of the roots of 
the trenscendentsl equstlons: 

sin s • ± * (D-la,b) 

end briefly describe s method of solution. 

Equations (D-la,b) both admit the trivial solution zero. They 
also have an infinite number of complex roots which are symmetrically 
lc'i’-red in each quadrant of the complex z plane. We will restrict 
our attention to the roots which lie in the quarter-plane (Re z > 0 , 

Im z < 0) , and all others will be obtained by symmetry with respect 

to the origin and the coordinate axes. Consequently, let us examine 
solutions of (D-la) of the form: 

z - a - i/3 a 0 ; j3>0 (D-2) 

Substitution of (D-2) into (D-ls) leads to the following system of 
coupled equations: 

sina cosh ■ a (D-3) 

cos CK sinhft * (S (D-4) 

Since cosh fi is necessarily positive for any value of )3 , (D-3) implies 
that sintt is also positive. Hence, possible values of a are such 
that 

2nrr £ a £ (2n+l)TT n • 0 , 1 , 2 , . . . (D-5) 

The above system may then equivalently be written in the following 
form: 
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F(a) = cos a a /( 



7 ^—) -1 - costi ( - /* -)= 0 
sina sinOT 


(D- 6 ) 


* ' C08h " l( slla ) 


(D-7) 


We have therefore reduced the problem to the determination of the 
zeros of F(a) . Such a function is only defined in the intervals given 
by (D-5). Its derivative is: 


dF _ (2a-sina ) 2 

da = 2 h 


(D- 8 ) 


Moreover : 

when a -» 2nrr+ and n^O F(a) -*+<*> 

when a 0+ F(a) -» 0 

when a -♦ (2n+l)n- F(a) — - <* 

It may be concluded from the resulting variations of F(a) sketched 

in Figures 19 that F(a) has only one zero CK + in each interval 

n 

(2nrr , (2n+l)rr ) # Since F( (2n+~)TT) < 0 each solution of (D- 6 ) can 
be further bounded as follows: 

2nn £ a* £ (2n+^)rr n ■ 0 , 1 , ... (D-9) 

n z 


and the corresponding value of j3 is then given by (D- 8 ) as: 


a 


= cosh *( 2 _) 


sina 


(D-10) 


With this information, one may immediately develop an iterative numer- 
ical scheme which will yield the roots 






(D —11) 


+ 

z 

n 




of the transcendental equation (D-la). The results are shown in Table 
III. Finally, as n goes to infinity in (D-9) and (D-10), both real 
and imaginary parts become infinite. From equations (D-3) and (D-4) 
we then obtain the following limiting form of z* : 

z + = (2n+|)TT - i cosh' 1 (2n+ ^)tt (D-12) 

n i l 

The transcendental equation (D-lb) may be studied in exactly the same 
manner* It is found that the roots are given by 

z" » a" - i/3" (D-13) 

n n n 


where 


(2n + l)rr < a ^ (2n +-)tt 
n l 


n = 0,1, . . 


and 


/3 = cosh (—7 — —=) 

n 'sma 

n 


•a 


(D-14) 


(D-15) 


The corresponding numerical results are displayed in Table IV. As the 
index n goes to infinity, the asymptotic limit of the roots is: 

z" = (2n +-|)tt - i cosh"*(2n+ r)rr 
n z l 


(D-16) 
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